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/ABSTRACT 

Seventy-five measurements of the differential elastic 

electron-deuteron scattering cross section were made for 

2 -2 
values of g ranging from 0.05 to 0.35 fm at the NPS LINAC 

with electron scattering angles of 60° to 120°. The charge 

radius of the deuteron was determined to be r„ = 1.9641 + 

E a 

0.0074 fm. Values of the neutron charge form factors were 

calculated using Lomon-Feshbach Models 1, 5, and 15 and 

relativistic corrections to the deuteron wave functions. 

Lomon-Feshbach Model 15 gave — ~ G £ (0) = 0.0200 + 0.0058 fm^, 

dq n 

in closest agreement with the value of ^ (0) determined 

dg n 

by thermal neutron work. A Rosenbluth plot of the data at 
2 

q =0.2 shows good agreement with theory. 
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I. INTRODUCTION 



A. OUTLINE OF THE EXPERIMENT AND INTRODUCTION OF FORM FACTORS 

The primary purpose of this worh was to measure the charge 

2 2 

form factor of the neutron, G (q ) , at low values of q (the 

E n 

square of the Lorentz invariant four-momentum transfer from 

the scattered electron to the neutron) in order to resolve 

differences among several earlier experiments [1,2, 3, 4] • 

Secondly, by comparing the G E obtained with the use of 

n 

several different deuteron wave functions to the G pre- 

n 

dieted by accurate thermal neutron experiments, one can select 
the S state and D state combinations of the deuteron wave 
functions that are in best agreement with experiment. ‘ 

Thirdly, measurements of the deuteron cross section at 
2 

low values of q permit the rms charge radius of the deuteron 
to be determined accurately. 

These three objectives were achieved by measuring the 
elastic electron - deuteron differential scattering cross 
section. The elastic electron - proton cross section was 
measured in order to normalize the measured deuteron cross 
section. The electron beam was produced by the Naval 
Postgraduate School linear accelerator (LINAC) . 

Since the primary purpose was the measurement of the charge 
form factor, it is useful to consider briefly here what a form 
factor is. In Section II a more detailed explanation will be 
given. The cross section,——, for a relativistic electron 
scattering from a nucleus that has zero spin and magnetic 
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moment and consists only of a point charge has been calcu- 
lated by Mott [5], The experimental cross section is re- 
lated to the so-called Mott cross section by 

(£ L P = ffi Ltt g2 ^ 2 ) 

where G is the form factor. If the nucleus has no structure, 

. 2 
l.e., it is a point, G = 1 at all values of q . If it has 

2 

structure, G = 1 for q equal to zero and G is less than one 

2 2 
for q greater than zero. A reason why G (q ) is called the 

form factor [3] is seen when the cross section for elastic 

scattering of an electron from a nucleus of charge Ze with 

finite extent is written [6] as 



(do ^ (do \ f f p (r) e iq * r d 3 r ~j 2 

Vd-n. /Exp Vd-n. /Mott Lj Ze J ' 



(1-2) 



where P(r) is the nuclear charge distribution and c^q is 
the three-momentum transferred to the nucleus. The term in 
brackets is the form factor G and is the Fourier transform 
of the charge distribution. If one calculates the theo- 
retical cross section for elastic scattering, it is seen 
that the charge distribution involves the product of the 
initial wave function of the scatterer at rest in the labo- 
ratory frame and the final wave function of the scatterer 
after it has absorbed the recoil momentum c^q. In Section 
II it will be shown that a relativistic correction to the 
deuteron wave function is required because of its recoil 
velocity. This is a small (a few percent) but important 
correction in this work. 
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If we consider elastic scattering of an electron from 
an infinitely massive (that is, a nucleus which does not 
recoil) nucleus of charge Ze Equation 1-2 will correctly 
describe the scattering; p( r) is then the static charge 
distribution. However, for light nuclei where recoil is im- 
portant, such as deuterium and hydrogen which were used as 
targets in this work, the intuitive concept of the form 
factor as the Fourier transform of the charge distribution 
is not unambiguous. For example, if the charge distribution 
of the proton were spherical in one frame of reference then 
in another frame the Lorentz contraction of the distribution 
would cause it to be nonspherical . To describe form factors 
in a way more useful in theoretical work, they are considered 

as invariant functions of the four-momentum q rather than the 

-* 222 22 ”*2 2 
three -momentum cftq, where ^ c q = ft c q - (E^-E^) and E^ 

and E^ are the initial and final energies. When the recoil 
2 ~*2 

is small then q ^ q and the concept of the form factor as 
the Fourier transform of the charge distribution has useful- 
ness. As will be seen in Section II the invariant form 
factors have other theoretical meanings in terms of strong 
interaction dynamics. 

There is another complication. Particles and nuclei may 
have magnetic moment distributions as well as charge distri- 
butions. These distributions may or may not have the same 
form as the charge distributions. Thus, a particle or nuc- 
leus that has a magnetic moment (but no quadrupole or higher 
moments) will have to be described by two independent form 
factors. The decomposition into electric and magnetic 
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distributions depends upon the frame of reference. For ex- 
ample, the elastic scattering of an electron from a massive 
nucleus with charge Ze and no magnetic moment will be cor- 
rectly described in the laboratory frame by Equation 1-2 if 
the nucleus is at rest in the laboratory frame. If, however, 
the calculation is performed in the electron's frame of refer- 
ence then the moving nucleus will also have a magnetic field 
and thus an apparent magnetic moment. In this work we are 
ultimately interested in the electric form factor of the 
neutron. A particularly useful frame, the Breit or "brickwall 
frame", separates the electromagnetic interaction of the nuc- 
leons into two form factors, G„, the electric form factor, 

ih 

and G^, the magnetic form factor. Sachs et al [7,8] showed 
that these are particularly useful form factors since G and 
G^ can be expressed as the Fourier transforms of the spatial 

distributions of charge and magnetization in the Breit frame. 

2 ~ } 2 

The Breit frame is the frame in which q = q [7], i.e., that 
has E^ = E f and that has the absolute value of the momentum 
of the incoming nucleon the same as the absolute value of the 
momentum of the outgoing nucleon and oppositely directed [9]. 

Consider the electric form factor of a nucleon in the 
Breit frame. 



which for a spherically symmetric particle may be written as 




o 
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2 

For small values of q the sine may be expanded to give: 

G E (q2) = G E (0) " 6 q2 r rms ' < I “ 3 ) 

where r rms is the root mean square radius of the charge 
distribution. A similar expansion may be made for the mag- 
netic distribution to obtain the magnetic r . The proton 
and neutron each have different form factors. Thus, four 

form factors are required to describe them: G e # G E ' G M ' 

P n p 

and . The ^ rms obtained by using Equation 1-3 is deter- 
n 

mined in the Breit frame, but since Equation 1-3 is only good 
2 

at low values of q , the r rms value determined is generally 
quoted as the nucleon's charge radius. 

B. METHOD OF OBTAINING THE NEUTRON ELECTRIC FORM FACTOR, G„ 

e n 

The process of obtaining G^, from the measured elastic 

n 

electron-deuteron cross section will be outlined here and 
discussed in detail in Section II. First, the elastic elec- 
tron-proton cross section was measured at a given scattering 
angle and scattered energy using hydrogen gas as a target. 

The hydrogen was then replaced with deuterium gas and the 
elastic electron-deuteron cross section was measured at the 
same scattering angle and with the same scattered energy. 

The incident electron energy was adjusted to give this con- 
dition. All geometries, beam currents, gas pressures, and 
operating procedures were kept as constant as possible while 
both the proton and the deuteron data were taken. In fact, 
the LINAC was operated continuously from the time an experi- 
ment started until it was completed with both cross sections 
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measured. The time required for such a "run" varied from 
about 15 to 40 hours. 

Because absolute cross sections cannot be measured here 
owing to systematic uncertainties, the ratio of the cross 
sections, cj^/qP, was formed. In this ratio, many systematic 
uncertainties cancel since the cross sections are measured 
under the same conditions. The absolute proton cross section 
has been measured in previous experiments [10,11,12], and 
these data have been fit by deVries [10] to formulas from 
which the cross section can be calculated. When the experi- 
mental ratio above is multiplied by the absolute proton cross 
section, the absolute deuteron cross section is obtained. 

Several steps are required to extract from this cross 

n 

section. The experimental cross section is first divided by 
the Mott cross section of the deuteron to obtain the deuteron 
form factor, G^. G^ accounts for all the structure in the 
deuteron, i.e., the electric and magnetic distributions 
within the neutron and proton as well as the convection of 
charge due to their motion within the deuteron and the mag- 
netic moment caused by this convection of charge. Following 
Jankus [13] , the deuteron was treated as though its two con- 
stituent particles, the neutron and proton, were point par- 
ticles. The distribution of charge predicted by the deuteron 

wave function is then contained in a term called the deuteron 

2 

charge structure factor, which depends only on q and the wave 
function of the deuteron ground state. E.L. Lomon [14] has 
provided three selected sets of values, from his larger group 
of sets of values, for these deuteron charge structure 
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factors. Each of these is based on a different deuteron wave 



function obtained from nucleon-nucleon data [15] . 

Drickey and Hand [3] experimentally verified a proportion- 
ality between the magnetic form factor of the deuteron and its 
charge form factor, called a "scaling law". That relationship 
was used to remove that part of the deuteron form factor 
caused by the intrinsic magnetic moments of the proton and the 
neutron as well as that part of the deuteron magnetic moment 
caused by the convection of charge. 

At this point only the extended charge distribution of 
the proton has not been accounted for in the extraction pro- 
cess. This distribution is accounted for by the proton elec- 
tric form factor derived from deVries' data (see Appendix A). 
The neutron electric form factor, G„ , has thus been extracted 

E 

n 

from the deuteron cross section. 

The above procedure is justified since the deuteron is a 
weakly bound structure, i.e., the proton and the neutron spend 
most of the time outside the range of nuclear forces and one 
can adequately describe the charge and current distribution 
within each nucleon as if they were free [13] . 

Early in this work, before the gas target system was 

operating, efforts were made to measure by a process simi- 

an 

lar to that described above but using polyethylene targets of 
CH^ and CD^ to obtain the cross sections. There were several 
problems with that work but the most important limitation was 
the scattering from carbon which complicated the data reduc- 
tion. The solid target data which were obtained were not as 
reproducible as the gas target data and are not included in 
this work. 
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C. RESULTS OE PREVIOUS EXPERIMENTS AND OE THIS WORK 

In this work the neutron charge form factor, G , was 

Jj 

2 n -2 
measured for values of q ranging from 0.05 to 0.35 fm 

2 

The slope of G versus q was obtained by a linear least 
E n 

squares fit to the data. 

2 

The slope in the limit that q 0 has been measured in 
many experiments [16,17,18] which measured the interaction of 
thermal neutrons with atomic electrons. The most accurate 
and most recent measurement of this kind is that of Krohn 
and Ringo [4] , who obtained 



-^2 G (0) = 0.0193 + 0.0004 fm 2 . (1-4) 

dq n 

This is probably that best determined of all nucleon form 

factor parameters [19] . Drickey and Hand [3] measured G„ 

E n 

by measuring the ratio of the elastic electron-deuteron 
cross section to the elastic electron-proton cross section 

using liquid hydrogen and deuterium for targets and found 

2 -2 
that G„ =0.00+0.01 over a range of q from 0.3 to 2.2 fm . 
E 

n 

That is, they found 



- A - g ^ 
, 2 

dq n 



0 , 



which is in disagreement with the thermal neutron measurements. 

Casper and Gross [20] developed relativistic corrections 
to electron-deuteron scattering and reanalyzed the data of 
Drickey and Hand. Using deuteron wave functions derived from 
the Feshbach-Lomon nucleon-nucleon model they found that the 
data of Drickey and Hand indicated 

G e (0) = 0.0153 + 0.0067 fm 2 , 
dq z n 
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which was consistent with the thermal neutron slope. 

Bumiller, et ajL. [21,2] extended the elastic clectron-deuteron 

2 2 -2 

scattering data to q of 0.10 ^ q ^ 0.8 fm , using CH 2# CD 2 , 

and C targets and substantiated the conclusions of Casper and 

Gross by obtaining a slope of 0.0179 + 0.0036. This quoted 

error of + 0.0036 is incorrect. It should have been quoted as 

+ 0.0089. Using the same equipment as Bumiller, Mader [1] 

continued the elastic electron-deuteron scattering work but 

used a gas target to avoid the carbon subtraction procedures 

which were necessary with the solid targets used by Bumiller 

et al . Although Mader ' s data are as accurate as Refs. 21 and 

2, his measurements showed that G ^ 0 over a range of 

dq n 

2 -2 

0.05 ^ q ^ 0.60 fm . In this work we have resolved this 

discrepancy by measuring the G- slope to be 

n 

-Ar G (0) = 0.0200 + 0.0058 fm 2 , 

0q n 

which then compares favorably to the value of 0.0193 + 0.0004 
2 

fm as measured directly by the neutron-electron scattering 
experiments . 
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II. THEORY 



A. ELASTIC ELECTRON-PROTON SCATTERING FROM A POINT PROTON 

Although it has been proved [22] that the proton is not a 
point particle, it is convenient initially to solve for the 
cross section for elastic electron-proton scattering from a 
point proton. 

The initial assumptions are: 1) the proton has a point 
charge and a point Dirac magnetic moment, 2) the one virtual 
photon exchange adequately describes the scattering process, 
and 3) the electron is a point particle. The Feynman diagram 
describing this process is 




where p^, p f and P^ and P f are the incident and final electron 
and proton four-momenta respectively and c^q is the four- 
momentum transferred in the scattering process, i.e., 

ft 2 c 2 q 2 = (p ± -p f ) 2 = {(Pj^-Pf) 2 - (E i -E f ) 2 }= (Pi-Pf) 2 . 

In order to obtain the cross section for this elastic 
scattering process we make use of S-matrix formulism. The 
S-matrix is a matrix connecting the initial and final states 
of the process shown above. The electromagnetic interaction 
responsible for electron scattering may be treated by per- 
turbation theory. The first approximation of the S-matrix 
is given by an integral over the interaction Hamiltonian, 
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S = 1 - i §dt H i ( t) , 

where H_^(t) may be written as a three dimensional integral 
over a Hamiltonian density H^(r,t) 

H. (t) = \d 3 r H. (r, t) . 

j. o i 

Si3ice the electrons can be described by a second quantized 
Dirac field and the proton's electromagnetic field may be 
described by a vector potential, the interaction Hamiltonian 
density can then be written as 

H ± (x) = - j^(x)A^(x) , 

where x is the four dimensional space time vector. The four- 
current of the electron can be written in terms of the four 
component state vector and Dirac matrices as 

J®(x) = - ie Y f (x) y^x), with ¥(x) = y + y 4 - 
The matrix element of this scattering, H f ^, is then 

H fi = - ie (x)y pi T i (x)A^(x)d 3 r (II-l) 

giving the interaction of the electron current (from its 
initial state to its final state) with the electromagnetic 
field of the proton. By making use of the first Born approxi- 
mation, the electromagnetic field due to the proton's transi- 
tion current from the initial to the final plane wave state 
can be obtained by solving the inhomogeneous wave equation 

(v 2 - ) A^x) = -4tt > 
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where the expectation value of the proton's current operator 
is the transition current creating the potential A. (x) . The 

r* 

solution is 

ATT — 

A^(x) = ~~2 ie 'Pf ^ x ) y M' Cp i ^ ' 

where <p (x) is the state vector of the proton. Substituting 
this result into Equation II-l the matrix element is 



H fi = - 4Tre 2 ^ f <x)y/.(x) V f (x)y Mi 9 i (x)d 3 x. 

q 



(n-2) 



Because of the use of initial and final state plane waves, 

the integration over space results in the delta- function 

^ (p^+P f - prp^) , expressing momentum conservation. It is 

2 

interesting to note that the 1/q term can be considered to 
be with either source current. Thus the proton may be con- 
sidered as being scattered from the electron’s field or the 
electron may be considered to be scattered from the field of 
the proton as it was above. 

The cross section may now be obtained by calculating 



da 

d-A 



U> 



E 

final 

states 



fi 

flux 



(H-3) 



and averaging over the initial spin states, where 



w fi = TT lH fi |6(E f - Ei ) 

and the flux is the incoming flux of electrons. The sum 
over final states is really an integral over all possible 
proton momenta and over that range of electron momenta which 
are within the solid angle of the detector. These integrations 
over momentum space and the averaging over the initial spins 
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and summing over the final spin states are tedious but then 
completed Equation II-3 yields 



= (I 

d .n/' \4 

Dirac 

proton 



e 4 cos 2 Q/2 \ / 



E? sin 4 9/2> 



1+ 2^' tan 



2 £ 

2 



)(- 



2E i sin k! 9/2> 



(II-4) 



1 + 



where E^ is the energy of the incident electron and 9 is the 
angle through which it is scattered. Both are measured in 
the system where the proton is initially at rest and 



c 2 fr 2 q 2 = 4E^E^sin 2 9/2 
E^ is the final electron energy given by 



E. 

E f 2E, ^ , ■ 

1 + — ~ sin 9/2 
me 

2 

In the limit that m^c /E approaches zero the first 
enclosed term on the right of Equation II-4 was calculated by 
Mott [5] in 1929, using the Dirac equation, to be the cross 
section of an electron scattered from a coulomb potential; 
in the case of a proton that physically means that the proton 
is a point charge and that it has zero spin and an infinite 
mass, i.e., it doesn't recoil. Let us call that the histori- 
cal Mott cross section. The final enclosed term on the right 
is the recoil term which accounts for both the kinematics of 
the recoil and the fields associated with the accelerating 
proton. When the historical Mott cross section is multiplied 
by this recoil term it is generally called the Mott cross 
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section. That nomenclature is followed in this work. The 



2 

tan 9/2 part of the second enclosed term accounts for the 



scattering contributed by the Dirac magnetic moment of the 
proton while the 1 in that term accounts for the point charge 
scattering. Equation II-4 thus provides the cross section 
for an electron scattering from a point proton with a point 
Dirac magnetic moment. However, there are two corrections 
which must be made to Equation II -4 in order that the theo- 
rectically predicted cross section of a proton agrees with 
measurements . 

The first thing which is known to be wrong with Equation 
II-4 is that the magnetic moment of the proton is not equal 
to one nuclear magneton as is predicted by the Dirac equation, 
but is measured to be ^ = + 2.79275 + 0.0003 nuclear mag- 
netons. Also the prediction of the Dirac equation of zero 
magnetic moment for the neutron compares badly with the 

measured neutron magnetic moment of ^ = -1.9135 + 0.0003 

V — 

nuclear magnetons. 

Pauli [23] has shown that the Dirac equation for a 
charged particle interacting with an electromagnetic field 
can be modified so as to represent a particle with an arbi- 
trary magnetic moment by addding the operator 



to the equation. Foldy [24] has verified that such a term 
satifies the restrictions of being Lorentz covariant and gauge 



form an infinite series involving arbitrarily high derivitives 




invariant and has shown that the general interaction terms 
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of the electromagnetic potential evaluated at the position 
of the particle. The series of coefficients of these terms 
can be interpreted as a series of moments of the charge and 
current distributions associated with the particle. We will 
return to this expansion later in this section. 

Because the four potential, A ,• can be regarded as the 
potential produced by the electron, the Pauli term can be 
shown to reduce to 



JC 



2m 



^2 (p i " p f ? 



which is added to the V term of the nuclear current. Here 
= 2.79- 1.0 = 1.79 is the anomalous magnetic moment of the 
proton in units of nuclear magnetons and 



O = _ Ify y _ y y ) 

M-V 2 H V v 



When the Pauli term is inserted into the matrix element of 
Equation II-2, the cross section for a point Dirac particle 
with its normal point magnetic moment and now also its 
anomalous point magnetic moment is obtained: 




point 

proton 



W/. 



\d svj l 1 + 

Mott 




2 2 , 0 . , 
tan (-j) + 




where 
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Mott 
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The second change which must be included in the description 
of scattering from a real proton will now be considered. In 
1955 Hofstadter and McAllister [22] measured the cross section 
of the proton and showed that the structure of a proton is 
more complicated than a point charge and point magnetic moment; 
the proton scattered as though it had a charge and magnetic 
moment "cloud" spread out to an rrns radius of about 0.7 fm. 
Subsequent work (1960-61) on the proton at Stanford [25,11] 
showed that the magnetic moment associated with the Pauli term 
had a softer or more spread out distribution with an rrns radius 
of about 1 fm, while the distribution of the charge had an rrns 
radius of about 0.8 fm. These rrns radii were calculated on 

the basis of Equation 1-3 by assigning proton form factors 
2 

Ffp (q ) to the Dirac charge and magnetic moment distribution 
2 

and F 2 p (q ) to the Pauli magnetic moment distribution. The 
addition of the, anomalous magnetic moment to the interaction 
and now the provision permitting the Dirac and Pauli components 
to have spatial extent changes the proton contribution in the 
matrix element of Equation II-2 from 



to 



v f (x)r^9i(x) 



'Pft.hp'^V WPiV-Pfv> F 2p^ 2 >)v 



The elastic electron-proton cross section with this 
modified current matrix element has been calculated by 
Rosenbluth [26] to be 

£ V2p) 2t » 2 <f>«p 2p y 

and is called the Rosenbluth cross section. 
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B. NEUTRON FORM FACTORS 

The neutron has zero total charge. Both the electron and 
neutron have magnetic moments, however, and thus they can 
interact through the magnetic dipole-dipole interaction. 

Also, an incident electron has a magnetic field associated 
with it because of its charge and this will interact with 
magnetic moment of the neutron. These interactions are well 
understood and must be accounted for in the calculation of 
the electron neutron scattering formula. 

There remains a third interaction which is generally called 
the neutron-electron interaction in the literature. This in- 
teraction is expected to be both spin and velocity independent 
if there is a separation of electrical charge in the neutron. 
Thus an electron, or any charged particle, penetrating the 
extended charge distribution of the neutron will be subjected 
to electrostatic forces. These may arise from two sources: 

1), the intrinsic charge separation in the neutron analogous 
to the spreading out of charge in the Dirac proton and evi- 
denced by the form factor F^ n analogous to F^ of the proton, 
and 2) , the anomalous magnetic moment of the neutron, des- 
cribed by a form factor ana l°U ous to ^2p ^ or P r °t° n / 

if the anomalous magnetic moment has a spatial extent. 

Foldy [27,28] has shown that since the neutron satisfies the 
Dirac equation, and since it possesses an anomalous magnetic 
moment, then some charge separation is to be expected. He 
compared the amount of charge separation predicted because of 
the neutron's anomalous magnetic moment with that charge 
separation inferred by the very accurate thermal 
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neutron-electron measurements of the neutron's attractive 
potential for the electron made by several groups [16,17,29,30]. 
The predicted and inferred charge separations agreed within 
experimental error, implying that the so-called intrinsic or 
Dirac separation is not present. Since this requirement for 
charge separation within the neutron simply because it has 
an anomalous magnetic moment is not intuitive, this will be 
examined in more detail below, as will the results of the 
experiments that measured the neutron's attractive potential. 

Several methods of measuring the attractive potential 
have been used. Only the first [16] of these is outlined 
here. If one measures the coherent scattering of neutrons 
by spinless (in order to remove effects due to magnetic in- 
teractions) atoms containing a nucleus and Z electrons, then 
the coherent scattering amplitude for neutrons whose wave- 
length is long compared to the size of the atom will be the 
algebraic sum of the scattering amplitudes due to the nucleus 
and the Z electrons. By decreasing the neutron wavelength 
until it is comparable to the size of the electronic cloud of 
the atom, destructive interference between the scattering from 
the individual electrons within the cloud of electrons can be 
achieved and the total coherent scattering cross section can 
be made to approach that of the nucleus alone. By separating 
these scattering effects, one obtains the so-called neutron- 
electron interaction. The result of the thermal neutron- 
electron experiments is expressed as the volume integral of 
the potential of the neutron-electron interaction. Plots of 
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the neutron (and the proton for comparison purposes) charge 
density and the interaction potential are given in Figure 
II-l [31] . 



p H NEUTRON 



Figure II-l 

Neutron and Proton Charge Densities and Interation 
Potentials Between These Particles and an Electron 

The neutron-electron potential is denoted as that 

is, it is the difference between what the neutron potential 
would be if it had no charge distribution and what it would 
be if it had a positive center and a negative outer region. 
This concept of a positive center- and a negative outer region 
is justified by the results of the thermal neutron-electron 
scattering experiments and is also consistent with meson 
theory. The volume integral of this charge distribution 
does not permit an independent determination of the depth 

and volume of the potential. It does however permit the de- 
2 

termination of the r of the charge distribution by the 
rms 3 2 

following calculation. 
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Let 



r 



.2 

rms 




and then use Poisson's equation to obtain 



r 



2 

rms 



^ ,^Vd 3 r. 



(XI-6) 



However, following Fermi, all experiments have reported 
results in the form of a potential V Q which was the average 
V which the electron sensed. This required that some aver- 



age radius of this potential be adopted. By convention, that 
radius is the classical electron radius (2.81794 fm) , which 
implies that classically the bound electrons in the atom 
were assumed to be interacting almost uniformly over the 
neutron, because as will be seen, the neutrons largest r 

rms 

is less than 1 fm. That is, the electrons are not just 
sensing the outer charge layer as might be implied by the 
classical picture of a point electron of low momentum pene- 
trating only the outer shell of the neutron charge distri- 
bution. This is perhaps best seen from the virtual photon 
scattering picture. The wavelength of the virtual photon is 
very long compared to the neutron diameter and the electro- 
magnetic fields will be effectively constant at any time over 
the extent of the neutron with the greatest difference in 
these fields occurring between the extreme edges of the 
neutron. Thus, while it is the charge distribution near the 
edges of the neutron which show the deviation of the neutron 
from a particle with zero charge density everywhere, it is the 
average charge density, and through Poisson's equation the 



o 
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average potential within the neutron, that is measured in 
these thermal neutron experiments. Using this average po- 
tential the volume integral of the interaction potential is 
then written 



S v <3 3 r = V o $ 



2 2 4 3 4 V o e 2 

4u r dr = ^ nrfv = - tt - 2 — r “ 
3 e o 3 2 e 

me 



(II-7 ) 



where r is the classical electron radius, 
e 

Since Equations II-6 and II-7 are simply different ways 
of performing the volume 'integral they must be equal, there- 
fore we have: 



thus 




me 



2 



r 



rms 




me me 



where V q is put in the conventional units of million electron 
volts. Then 

r rms = 31 * 078 V o (fm)2 • (II-8) 

The presently accepted value of V q is [4] : 

V = 3720 + 90 electron volts . (II-9) 

o — 

This value of V q when substituted into Equation II -8 gives 

2 2 
r rms c ' har 9 e of neutron = - (0.34 fm) , the negative sign 

showing, since the charge on the outer shell is weighted by 
2 

r more heavily than the inner shell, that there is a charge 
separation in the neutron, and that the outer regions of the 
neutron are negative. 
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It was stated above that the two facts that the neutron 
satisfies the Dirac equation and that it has an anomalous 
magnetic moment are sufficient to account completely for the 
measured value of the potential V o [27,31,32] . Let us now see 
why that statement is correct. 

The relativistic Hamiltonian for a neutron with the 
anomalous magnetic moment as attributed by Pauli can be 
written as 

H = 0M + S*P-3f n (e/2m) -.i/sS-Ej 

where K n is the anomalous (and total) magnetic moment of the 
neutron measured in nuclear magnetons. This can be reduced 
to the corresponding nonrelativistic Hamiltonian by the method 
of Foldy and Wouthuysen [33] which yields 

H = AM + 0p 2 /2M - K div E 

n 4tT 

+ K [P x E - E x P] 

n 4M 

2 

with terms of order greater than (1/M) being neglected. 

For the coulomb field of a point electron located at 
the above Hamiltonian becomes 

H = $M+ (/3p 2 /2M) + 4 tt 3C ( e 2 /3/4M 2 ) 6 (x-x ) + . . . 

where the term containing the delta function is in the form 
of the electron-neutron interaction. Equating the volume 
integral of the neutron-electron interaction potential to 
this term one has [27] 

V o = ^\(e 2 /M 2 ) [ (4TT/3) (r 2 )] -1 = - 3900 eV. 
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This compares very well with the experimental result given in 
Equation II-9. Since the neutron charge distribution can be 
completely accounted for by its anomalous magnetic moment this 
means that there is, within experimental errors, no Dirac 
particle intrinsic charge distribution unless relativistic 
contributions fortuitously cancel this Dirac intrinsic charge 
distribution. 

The distribution of the neutron anomalous magnetic moment 

was measured by Yearian and Hofstadter and Bumiller [34,35] 

in 1958. They used Rosenbluth's proton cross section as given 

in Equation II-5 with modifications to obtain the neutron 

cross section. Based principally on the thermal neutron work 

discussed above, they set F^ n = 0 to describe the neutron's 

lack of intrinsic charge and magnetic moment distribution and 

allowed to describe its anomalous moment. The neutron 

zn 

cross section was then 



dcr = AdcrN 

d-n. \d^v/ 



Mott 



2n n 



4M 



Electrons were scattered from a deuterium gas target at suffi- 
ciently high energies and large angles to enhance the magnetic 
scattering from the nucleons. By using the rather good approx- 
imation for inelastic scattering from a deuteron, i.e., 

IN 

tf = a 4 - cr they were able to obtain the r of the ano- 
d p r\ J rms 

malous magnetic moment of the neutron and concluded that it 
was between 0.8 and 0.9 Fermi. 
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C. CHARGE , MAGNETIC AND ISOTOPIC FORM FACTORS 

Thus far the four form factors used to describe the dis- 
tributions of charge and magnetic moment of the Dirac particle 
and the Pauli moment have been introduced in the historical 
manner in which they evolved. They may be summarized as 
follows : 

F n is the form factor associated with the distribution 
IP 

of charge of the Dirac proton and also with the distribution 

of the intrinsic (Dirac) magnetic moment. In the static 

(q 2 = 0) limit F (0) = 1. 

is the form factor associated with the anomalous 
2p 

(Pauli) magnetic moment distribution of the proton. In the 
static limit F 0 (0) = K = 1 . 79 . 

2p p 

F n is the form factor associated with the distribution of 
In 

charge of the Dirac neutron and also with the distribution of 
the intrinsic (Dirac) magnetic moment of the neutron. In the 
static limit (0) = 0. 

is the form factor associated with the anomalous 
2n 

(Pauli) magnetic moment distribution of the neutron. In the 

static limit F 0 (0) = 3C = 1.91. 

2n n 

It has been shown above that the anomalous (Pauli) mag- 
netic moment of the neutron causes some charge distribution 
within the neutron. Therefore the form factor not only 

tells what the anomalous magnetic moment distribution is, but 
tells something about the charge distribution. 

The work of Yennie, et al. (1957 [31], Ernst et al. (1960) 
[8] and Hand et aT. (1963) [19] introduced linear combinations 
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of the four Dirac and Pauli form factors for both theoretical 
and experimental reasons. These combinations were the 



following: 



2 



G. 



E 



P 



F 






In Ref. 19 it is pointed out that these linear combinations 
correspond to zero, G E , and one, G^, unit of angular momentum 
transferred along the direction of the virtual photon ex- 
changed in the scattering process. Reference [8] added more 
physical content to these linear combinations of and F£ by 
showing that the rms radii of and charge and magnetic moment 
spatial distribution are derived from Equation 1-3 by 



G^ are the "charge" and "magnetic" form factors. Additionally, 
Ref. 19 points out that the cross terms of the charge and mag- 
netic form factors do not interfere with each other in the 
Rosenbluth cross section when it is written in terms of G„ 

hj 

and G.., while the Dirac and Pauli terms are mixed into the 
M 

Rosenbluth cross section as products of each other, as can 
been seen by examining Formula II-5. Written in terms of G £ 
and G^, the Rosenbluth cross section for the proton becomes 




and 




at q‘ 



.2 



0. For this reason the descriptive names for G £ and 



34 



1 



acr 

d-TL 



m [- 



_Ep_ 



Mott 



1 + 



4m 



+ 2 tan 



0 / 2 ) G M J 



(11-10) 



i+ 



4m 



4m 



with an analogous equation for the neutron cross section. 

G e and are obtained experimentally by measuring the cross 

section at at least two different scattering angles, at the same 

2 2 
q . From the intercept of a plot of G versus the coefficient 

2 

of G fc/ , one obtains G„ and from the slope one obtains G.„ . 
ME * M 

P 2 p P 

The static (q =0) values of these charge and magnetic form 



factors are: 



G_. (q 2 = 0) = 1 



G m (q 2 = 0) = 2.79 (II-ll) 

P 

G e (q 2 = 0) = 0 



G m (q 2 = 0) = - 1.91. 



We now consider what the slope of G £ tells about the 

n 

charge distribution of the neutron by studying the primary 
scattering event of interest in this work, the e-d scattering. 



This process may be described by the Feynman diagram 
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where and are the incident and scattered electron 
states, d^ and d f are the target and recoiling deuteron 
states and q is the four-momentum transfer. The deuteron has 
a diameter of about 4 fm; in this work, with electron energies 
of about 35 to 100 MeV, the virtual photons had wavelengths 
in the range of 10 to 30 fm. With these long wavelengths 
compared to the extent of the deuteron and neutron, only the 
outer charge distributions are measured. It is important to 
note that the thermal neutron experiments [18,4,16,17] give a 

positive slope for --- g e as do two of the e-d experiments 
dq n 12 

[20,21]. This means that the quantity - -g- ^ rms is positive, 
i.e., that the outer part of the neutron is negative. This 
is so because the square of the rms radius is obtained from 
the expression 

CO 

r rms = Ji r 2 (4Trr 2 p(r)dr) 
o 

c° 

where the total charge of the neutron = \ 4rrr p(r)dr = 0; 

2 ° 

thus a negative ^ rms shows that the charge density, p(r), 
of the neutron must be more negative at larger radii. 
Historically, theoretical attempts have been made to explain 
nucleon charge distributions in terms of a third set of nucleon 
form factors which are called the isotopic form factors. These 
isotopic form factors are useful because the structure of nu- 
cleons is determined by the strong interaction and calculations 
can be arranged in such a way that one can take advantage of 
the conservation of isospin. As the form of the nuclear elec- 
tric current operator [36] 

1+ T 3 

Jp, (x) = e<p(x)yp, — 2 ‘P(x) 
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shows, the current has mixed transformation properties under 
rotations in isospin space; part of it transforms as a scalar 
and part as a third component of an isovector. As a conse- 
quence of these mixed transformation properties, the proton 
and neutron current elements may be expressed in terms of an 
isovector and isoscalar current, i.e., 

J*(x) = j®(x) + <(x) and 

J “(*) = J j>) - 

The electric and magnetic isotopic form factors then 
enter_ the cross section calculation as the coefficients of 
the scalar and vector current operators in the matrix element. 
The isotopic form factors, where the superscript denotes 
scalar or vector, are: 

4 = * [ g e + g e ] h 

p n ^ 

G E = ^ £ G E G E ] q 2 -* 0 ^ 

p n ^ 

(II- 12 ) 

Gm = % ] ~2~r > h (M- +M- ) 

H ^ t M MJq^O ^ ' p n 

p n ^ * 

k “ G. a "| 2 % (M* -M- ) 

M ^Lm M J q 0 ^ v p n 

p n ^ r 

There have been attempts to give the isotopic form factors 
further meaning. While these attempts have not been completely 
successful, it is useful to consider them briefly. A state of 
n tt mesons, produced by a photon, must have zero total charge 
and be odd under charge conjugation. Based on this, Drell and 
Zachariasen [37] have shown that only an even number of pions 
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contribute to the isotopic vector form factor and that only 
an odd number of pions can contribute to the isotopic scalar 
form factors. They also have shown that charge conjugation 
does not permit a single pion to be produced since the photon 
is odd under charge conjugation and the tt°, the only tt which 
conserves charge, is even under charge conjugation. The 
possible intermediate states then include: 





Figure II-2 

Possible Intermediate Photon-Proton States 



In the diagram in Figure II-2 the possible existence of 
pion-pion forces were disregarded; but we know that there are 
resonances in the two-pion and three-pion systems, and we 
expect that such resonances may have an effect on the elec- 
tromagnetic structure of the nucleons. Chew et_ al . [38] and 

Nambu [39] have shown that such propagators as depicted in 
Figure II-2 (b) and (c) do in fact contribute to the form 
factors in such a way that they have poles at the masses of 
the P q , and The rho (isospin = 1, spin = 1) contributes 

to the isovector and the omega and phi (isospin = 0, spin = 1) 
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contribute to the isoscalar form factors. On the assumption 
that virtual clouds of mutually strongly interacting bound 
pions are responsible for the electromagnetic properties, 
these clouds show up as resonances in dispersion theoretical 
expressions for the form factors, for instance 



<^(q 2 ) 



_ 03 O P 

l C 

n J_ _ . 2 2 



(II-13) 



where m corresponds to the threshold for resonance creation. 

Furthermore, if one assumes that these vector mesons (spin =1) 

2 

appear as sharp peaks in the spectra of weight functions g(q' ), 
i.e., neglecting the widths (P (100 MeV) , to(10 MeV) , <p(3 MeV) , 

Equation 11-13 becomes 



G®(q 2 ) = |[(l-S el )+ ], (H-14) 

1 + 

M S 

which is the form factor expression which Clementel and 
Villi [40] have proposed, where is a constant. The term 

1-S 0 ^ represents the possibility of high mass resonances, 
i.e., it represents those resonances which have such high 
masses that their poles may be replaced by constants. M g is 
the mass of the scalar resonances. There is one term for 
each resonance. Similar expressions for G^, and were 
likewise developed and were used by deVries [10] to fit proton 
and neutron data. This fit is discussed in Appendix A and 
provided the normalization for the proton data obtained in 
this work. 
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Earlier in Section II it was observed that the outer 
regions of the neutron have been measured to be negative in 
the thermal neutron work. There are several theoretical 
reasons to believe that the outer regions should be negative 
[32, 38], Among these is. the fact that the anomalous magnetic 
moment is almost entirely a vector in isotopic spin space, 
i.e., the anomalous moments of the neutron and proton are 
nearly equal in magnitude with opposite signs. Additionally 
since the isovector meson, the rho, is the least massive of 
the rho, omega, phi mesons, the uncertainty principle as well 
as the Yukawa potential indicate that the outer regions will 
be dominated by the rho meson. 



D. METHOD OF EXTRACTING THE PROTON ELECTRIC FORM FACTOR FROM 
THE PROTON CROSS SECTION 

The Rosenbluth cross section for the proton (Eq. 11-10) 
can be written as: 



(d£) = (d£) 

vw . W/. 

proton 



[iTr + if? g m + 2Tg m tan " 

Mott p p 



!]• 



(H-15) 



2 2 

where T = q /4M^. This can be simplified by using the scaling 
law proposed by Lehmann et ad. [41] and others [19] , which 



states that 



M„ 



(ll-15a) 



Schumacher has shown that we should not expect isoscalar 

form factors to equal isovector form factors [42] , and that a 

set of scaling laws which relate the form factors and G*, 

MM 
p n 



to G,. 



and G,_ 



G„ 

P E 



^n G E 
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and 



M 



M* G + M- G . 
n E p E 
p n 



For our experiments, magnetic scattering generally accounted 
for only a percent or two of the total scattering in most runs 
and at most 10% of the scattering in a few runs; thus the dif- 
ference in these two scaling laws, which was at most a 1.0% 

2 

difference at 0.35q , then made a difference in the cross 
section of less than a nominal 0.01% ux^ to at most a 0.05%. 
Thus the conventional scaling law was used to reduce the cal- 
culations. When the magnetic contribution is accounted for 
by using the scaling relationship. Equation II-15a, Equation 

11-15 reduces to 
^2 



dof 

d-n. 



E p 

°m (0) ITr L 1+ 7|J 'p + ( 1 + T ) (2T)n2 tan 2 |J < 



(11-16) 



where 



acj p 

d.n 



Mott 



= a p 
M 



(9). 



Because the last two terms in the bracket are due to the 
magnetic form factors it is convenient to group them together 
by defining 

C p = rv?+ 2T>i. p (l+T)tan 2 § . (11-17) 

Because the proton cross section can be written as 



&21 

d-n. 



= *>)G 2 , 



(H-18) 



where is the proton form factor, we can obtain the proton 
electric form factor by combining 11-16, 11-17, and 11-18 as 




.2 (1+T) 



p(i+c p ) 





da p 

da. 

a M< S > 



(11-19) 
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where 



K 



P 

M 



u+n 

( 1+ V 



is defined as the proton magnetic correction term. Thus by 
measuring the proton cross section and dividing it by ^ (9) 
and then multiplying this quotient by one extracts the 
square of the proton electric form factor. 



E. THE DEUTERON GROUND STATE 

Before the elastic electron-deuteron calculation is studied 
it is useful to consider the ground state of the deuteron. 

This is particularly true in this case of elastic scattering 
because the deuteron is in its ground state both before and 
after the scattering event. It is an experimental fact that 
the deuteron has but one bound state and that is the ground 
state. Since the discovery of the deuteron in 1932 many of 
its characteristics have been measured and provide useful 
guidance in interpreting the electron deuteron scattering data. 
Some of the more important are discussed below. 

1) The neutron and proton are bound together by an energy 
of 2.2256 + 0.0015 MeV. 

2) The angular momentum is 1^, and the parity is even. 

3) The magnetic dipole moment is = 0.857393+0.00003 
nuclear magnetons. 

4) The electric quadrupole moment is 0.00282 + 0.00002. 

5) The root-mean-square electromagnetic radius of the 
deuteron is 1.96 + 0.02F. 

6) The deuteron potential has a repulsive core with an rms 
radius of about 0.4F. 
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7) The mass of the deuteron is 1875.585 MeV. 

Elementary calculations solve for the deuteron wave 

function in a square well neutron proton potential using the 

measured binding energy. The solution shows that the neutron 

and proton spend most of the time outside of the square well 

potential. The ground state is mostly a IS state, i.e., an 

l - 0 state, which is indicated not only from general energy 

considerations [43] , but from the fact that the algebraic sum 

of the proton and neutron magnetic moments nearly equals the 

deuteron magnetic moment. This indicates that the intrinsic 

spins of these particles are aligned and that there is no 

orbital angular momentum. If the nuclear force were spin 

independent, one should expect to observe a singlet state, 

1 with the same energy. However, no bound state with J = 0 
b 0 

has been found in the deuteron. This indicates that the 
neutron-proton force is stronger when the spins are parallel 
than when they are antiparallel; however, this does not mean 
that the antiparallel or singlet force is zero. In fact, the 
singlet force is almost strong enough to produce a bound state. 

Although the above description of the deuteron is accurate 
to first approximation, a better description is possible. 

First, one notes that the algebraic sum of the proton and 
neutron magnetic moment is 0.87925 + 0.0004 while the measured 
magnetic moment of the deuteron is 0.857393 + 0.00003 nm. 

This difference, about 2^f/ Q , is larger than the experimental 
error. The second problem with the description given above is 
that a pure S state implies that the quadrupole moment is zero 
since the wave function is spherically symmetric, while the 
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deuteron in fact has a positive quadrupole moment. The only- 
other state in which two particles give the same angular mo- 
mentum and parity is the " D-^ state = 2, S = 1, J = 1) . 

The conclusion is that the total wave function of the deuteron 
must have the form Y = a^Y^ + where the subscripts indi- 

cate S state and D state and the coefficients satisfy 
2 2 

a g + a D = 1 in order to normalize the total wave function. 

The exact value of a^ that produces a quadrupole moment in 
agreement with the experimental value depends on the radial 
part of the function Y # which in turn depends upon the de- 
tails of the radial dependence of the neutron-proton potential. 

2 

Values of the D state probability, = a^, in the range of 

3 to 1 % have been obtained from calculations using different 

2 

potentials [44] . A mixed deuteron ground state with a^ ^ 4% 
produces an electric quadrupole moment and a magnetic dipole 
moment, both in qualitative agreement with observations. It 
is one of the objects of this work to help determine which 
wave function best fits these values and from this determine 
more properties of the nuclear potential. Different wave 
functions which are considered in this work have different 
percentages of D state. However it is possible to keep the 
same percentage of D state and to change the quadrupole moment 
by changing the range of the tensor potential. This potential 
cannot be spherical as that produces a pure S state. The 
positive quadrupole moment indicates that the two nucleons 
prefer to line up with their spin vectors along the same line 
rather than side by side. The tensor potential which is a 
function of the angle between the spin vector of the neutron 
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and the proton and the radius vector r separating them accounts 
for this preferential alignment of spins and is a potential 
which must be added to the dominant spherical potential. 

Several textbooks [43,45,46] cover in detail all the con- 
cepts that were briefly discussed above. Having summarized 
the deuteron ground state we now consider in detail the 
elastic electron deuteron scattering process. 

F . ELASTIC ELECTRON-DEUTERON SCATTERING 

Jankus [13] in 195$ derived the cross section for elastic 
electron-deuteron scattering in the first Born approximation. 

In his calculation the nucleons were represented by point 
charges and point magnetic moments. The electron was con- 
sidered extremely relativisitic, but the deuteron was treated 
nonrelativi stic ally. 

While the nucleons have spin one half and two form factors, 
the deuteron, with spin 1, has three form factors. This is 
consistent with the work of Glaser and Jaksic [47] and others 
[48] which predicts that the cross section for scattering a 
relativistic electron from a potential of spin J will contain 
2J + 1 form factors. 

Jankus' s result, written in the more current terms [49], 
and including form factors for the proton and neutron, i.e., 
permitting Jankus' s point nucleons to have extent, is 

= ff M (0) { G o (q2)+ | r ) 2 G2(q 2 )+ |l7G^(q 2 )[l+2 (I+ 77 ) tan 2 f]}, (11-20) 

where TJ=q 2 /4M d and Gq (q 2 ) , G^q 2 ) and G 2 (q 2 ) are the charge, 
quadrupole and magnetic moment form factors of the deuteron. 
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These are given by 



G o = (G E. + g e )C E' 
P n 



(11-21) 



where 



C E = § [u 2 (r)+ w 2 (r)]j o ( q f)dr (11-22) 



for the charge from factor; 



G 2 = < G E + G E )C Q' 



p n 



where 



Cq = ~ \ [u (r) w (r) - -- 

Q WV J 0 L 2>I 2 



]j 2 ( q 2 ) 



dr 



for the quadrupole form factor, and 



(11-23) 



(11-24) 



M d 



{( G M + g m ) c s + ( G E + g e ) C J 



G 1 = 2 M~ U G M + G M ^ C S + V G E + G E 
P P n P n 

for the magnetic form factor, where 

. tK ( r )-!- 2 

o 

and 



(H-25) 






C L ■ I S “ 2 < I )[i 0 ('! |)+ j 2 (q f)]ar- (11-27) 



(G*„ + )C„ determine the contribution of the intrinsic 

M ■ M S 
P n 

magnetic moments of the proton and the neutron to the cross 
section and (G^ + G„ )C_ determines the magnetic contribution 

£» Jj Ju 



P n 

of the convection of charge in the deuteron to the cross 
section* 

Here M and M, are the mass of the proton and deuteron 

P Q 

respectively and u (r) is the S-state wave function and w(r) 
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is the D-state wave function normalized so that 

§ [u 2 (r)+ W 2 (r)]dr = 1 (11-28) 

o 

and j 0 (q ~) and j 2 (q -j) are spherical Bessel functions. 

The static limits of the functions above are 

s G q (0) = 1, G 2 (0) = M 2 Q, 

G l (0) = iT [V V f W V *>]' 

P 

h[G (0)+ G (0)] = G (0)= h.hiG (0)+G (0))=G m (0)=% (»+!*) . 

p n S p n S ^ 

C E (0)= 1, C Q (0)= M 2 Q, C s (0)= 1 - | P d (XI-29) 

and o 

C L^°^ = ~2 P D # 

Q is the quadrupole moment of the deuteron, M* and M* are 

p n 

the proton and neutron magnetic moments, and 

CO 

P D = J ^ ^ r ^ dr 

is the D-state probability of the deuteron. 

Introducing the nucleon form factors to Jankus's results 
of course assumes that the nucleons are not distorted in the 
deuteron. The fact that these nucleons are so weakly bound 
that they spend a sizable portion of their time far apart 
makes this assumption of non-distortion plausible. 

Gross [50] and Casper [20] have calculated relativistic 
corrections to the electron deuteron scattering problem. 

While these corrections could come primarily from relativistic 
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modifications of (a) the deuteron wave functions and (b) the 
nucleon current, the latter have been shown to be small by 
Gourdin [71] and also by Gross. The relativistic modifications 
of the deuteron wave function depend not just on the Lorentz- 
contr action factor, which is actually a contraction in two 
directions [31], but also on the fact that the equal-time wave 
function in the rest frame is not the same as the equal-time 
wave function in the moving frame [20] . Gross goes on to 
solve for the corrections to the deuteron wave function and 
then obtains workable expressions for the corrections in terms 
of u (r) and w(r). Gross's expression for G^, G^ , and G 2 are 
then much more complicated than those of Equation 11-21, 11-23, 
and 11-25. For example. Gross's complete result for the deu- 
teron charge form factor is 

C E =( 1 “ q 2 /32M 2 )^ (u 2 +w 2 )j Q (q “)dr- { (u 2 -uru' 

o 24M o 



- uku + 7w 2 - w r w’ - wk w ) [ j 0 (q[ f-)+ j 2 (<2 f")] dr 
2 ~ 

+ 2 \ (uru'+ wru' 4- wrw' ) j 0 (q|-)dr, (11-30) 

16M, °o 
d 



where prime refers to differentiation with respect to r and 

2 ["^2 -I 

k = r — j " a = l Me / and e = deuteron binding energy. 



The first integral is the usual Jankus nonrelativistic result. 
Schumacher and Bethe [44] and Elias et ad. [49] note that 

the relativistic corrections of Casper and Gross [20] to the 

2 2 

deuteron form factors, to order q /2M , enter simply as a 

P 

2 

factor (1- T/2) on the right hand side of Equation 11-20, 
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? 2 

where T - q"/4M . Equation 11-20 can be further simplified 
2 

for small q by using the scaling law [42] 



G i = TT 

1 d 0 M 



(II-31) 



which is not a law but an approximate fit to the data which 

can be derived theoretically from the nonrelativistic quark 

model and SU6 . (SU6 also predicts /M* =1.5 which is not 

p n 

quite true.) 

By making the relativistic correction as noted by 
Schumacher above and using the scaling law. Equation 11-20 
becomes 






( G o (qZ 



8 2 2 , 



2,,2 2 



)+ eft G 2 (q )+ 3 ' M, d G o T L 1+2(1+r?)tan 



2 £ 

2 



]}u- }) 



(11-32) 



Rewriting this using Equation 11-21 and dropping the quad- 

-4 

rupole term as it is less than 5 x 10 at the highest value 
2 

of q reached in this experiment, we then have 

= <^(0) (G +G E ) 2 C^{l+ §M. 2 r (1+2 (1+77) tan 2 |)}(1- j) 2 . (11-33) 

p n 

If we define C d = -|t|J. 2 (1+2 (1+7? tan 2 -|) ) (11-34) 

and note that 



then 9 

Additionally, 



= a^(9)G 2 , 

d-n. M d 



(11-35) 



2 2 , 



(G E + G E } S^d^ 1 - 2> 



T\ 2 



(11-36) 



if we define as the deuteron magnetic and 
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(11-37) 



relativistic correction we have: 
K mr = 'TTcT) ( ~r 

Cl 1- TT 



We obtain the so-called deuteron electric form factor G £ 
from Equation 11-36, i.e.. 



g e = (g e +g e )2c e 

d p n 



Gd j d 

’ (1+C d ) (1- f) 2 G d K mr~ K mr (0) 



. (11-38) 



With both the deuteron and proton electric form factors avail- 
able we now obtain the neutron electric form factor by first 

forming the ratio of G to G , i.e., 

a E p 



( g e (Od 5 + g ^ ( ^ ))c 

d _ P 



g e (< ^ } 
p ^ 



V’ 3 E 



G E * g p^ 
P F 



(11-39) 



2 2 

where the q, and q indicate that the proton and neutron m 

d p 

the deuteron have received a four-momentum transfer different 
from that of the proton which was the nucleus of the hydrogen 
target. This was so in this experiment because the incident 
electron energy was such that the recoil electron from both 
hydrogen and deuterium targets had the same final energy. 

Thus the transmission efficiency of the spectrometer and the 
detector efficiencies were held constant for both targets, but 
the four-momentum transfer was slightly different. Since the 

proton electric form factor is available from deVries [10] b' 

. 2 

fit we may scale the G £ (a ) measured from hydrogen to the 
2 p ^ 

G E ( q d^ measu red in deuterium. Thus the electric form factor 
P 
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of the proton in the hydrogen that would be predicted from 
the hydrogen measurement, if the hydrogen had experienced the 
same four-momentum transfer as did the proton in the deuterium 
is given by 



3 e ^ q a^ 
p 



g e < q p> 
p ^ 



3 e (q a^ 

p 

g e (q p> 
p ^ 



deVries 



(11-40) 



Equation 11-39 may be rewritten in terms of the measured 
quantities, the deuteron and proton cross sections, by using 
Equations 11-19 and II-3&. This gives 



[ g e (q a )+ g e (q a> 

L p n - 



G E ^ q p^ 
P 




(11-41) 



Now forming this ratio, i.e., Equation 11-41, with the 
same four-momentum transfer for the proton in both the 
deuterium and hydrogen we have: 



[ g e (q a )+G E (q a>] c e [f-- aP( 0) K d Y 

L p n d d-a m mr 



g e < q a> 



which gives 

g e < q a> 



( 



a d (8l kP 

a m m 




(11-42) 



deVries 



g e (q a> 

p 



<J p (8) k' 

m mr 

cr d (9) K P 
n m 




1, (11-43) 



deVries 



'E 



where R is the experimental ratio of the deuteron cross section 
to the proton cross section. 
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Note that cr p (9) and K p were calculated at the incident 
m m 

energy of the electron on the hydrogen, which was higher 

than the energy of the electron on the deuteron, at which 

0^(0) and were calculated, 

m mr 

2 

The neutron electric form factor, G^, (q^) / then is 

n 



given by 






<J P (Q) K° 
m v ' mr 

0- d (9) K P 
m m 



V 4 ’ 



deVries 



c - 1 *° 

C E 



V 4 ’ 



(11-44) 



which may be rewritten as 



1 7^" 



E n d 



G E C E 



1.0 



G E (q d^ ' 
P 



(II-44b) 



G q, is the G„ which includes Gross's relativistic 
E d E 

n n 

corrections. It appears at first to be quite different from 
those of Stewart and Mader because they used the nonrelati- 
vistic deuteron wave function to calculate G„ and then added 

E 

n 

the change in G„ which is due to this relativistic correction 
E 

n 

in the deuteron wave function to the nonrelativistically cal- 



V 

which is practically model independent within the range of q 
considered here [2], i.e.. 



) 



= G„ + Ag 



g e _+ 



iJ n rel n E n - n 
This can be shown to be very close to the correction of 
Schumacher used in this work by considering the following. 
Schumacher's relativistic correction enters as a factor 



(1 - T/2) on the right hand side of Equation 11-20. This 
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correction is carried forward to Equation 11-43. If we 
form the nonrelativistic ratio of G to G we have 




K d = K d 
mr m 



(1- T/ 2) 4 



- 1 . 



(11-45) 



If K were inserted in Equation 11-44 in place of K we 
mr m 

would have 




and 

AG. 



E 




8m 



(11-46) 



(11-47) 



(11-48) 



2 2 

where T = q /4m has been written explicitly. Since the 
2 2 

coefficient of q /8m^ in Equation 11-48 is experimentally 
within 2 percent of unity throughout this experiment, then we 
see that Schumacher* s method and Stewart's and Mader's method 
are in good agreement. 
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G. THE ROSENBLUTH PLOT 

Rosenbluth derived the cross section for electron-proton 
scattering based on the one photon exchange. Jankus made the 
same assumption in calculating the deuteron cross section. 
Examining Jankus's result, given in Equation 11-20, and, re- 
produced below. 



do 

d-a 



/Mott 



G 2 (qV f^G^qV frjG 2 (q 2 ) [ 1+2 (1+T?) tan 2 f] 



which can be written as 



do 

d.a 



do 

d-rt 



Mott 



(II-48a) 



(II-48b) 



it can be seen that if the deuteron cross section is measured 

2 2 
at different angles but at a constant q , then a plot of G^ 

2 0 

versus (1+2 (l+?7) tan -~) should give a straight line with 

2 8 2 2 2 2 
Gq + -g?? as the intercept and -^G^ as the slope. Thus a 

linear plot of such data would indicate that the one photon 

process accurately describes the phenomena. Such a plot has 

been made in Section VB. 

Equation II-48b cannot be used directly because that 
equation is presumed to give an absolute cross section. The 
absolute deuteron cross section is obtained through knowledge 
of the hydrogen cross section. Since deVries b' fit is the 
best fit to the proton absolute cross section, this is used to 
normalize the data. That is 

(^Absolute" (^deVries" " ^ ) LINAC (II_49) 

where N is the normalization factor for a given point. 
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Therefore 



N 



- ^ 
'. a * l 



da p 

dji 



(11-50) 



/deVries / \ / NPS LINAC 

Because the deuteron and proton cross sections were measured 
under the same experimental conditions, the absolute deuteron 
cross section is 



dcr 

d-a 



= N 



Absolute 



da 

d-n- 



NPS LINAC 



It follows that 



\ d - h 



Absolute 



doP 

d-n. 



deVries 




LINAC 



In terms of the experimental ratio R, where 



dcr" 

d-n- 



NPS LINAC , 



dcr p 

d-n. 



NPS LINAC 



da 

d-n. j 



But 



Absolute 



-) 

/Absolute 



daf 

d st 



da" 

dsL 



deVries 



Mott 



<■ 



therefore 

J2. 



where 



(0) 



da* 

d-n. 



deVries 



deVries 
and 



da 

d-n. 



R / 
m' 



(H-51) 






Mott 



MOtt y 

were calculated by programs 



T0PR0T0N and TOPDEUTR which are listed in Appendix E and the 
ratio R m was the weighted mean given by Equation IV-20. 
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The statistical error in G, was calculated from 

d 




where is calculated from Equation IV-21. 

R m 

Several points were taken in order that the data could 
be checked for consistency on a Rosenbluth plot. The results 
are given in Section V-B. 
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III. EXPERIMENTAL EQUIPMENT AND PROCEDURES 



A. THE NAVAL POSTGRADUATE SCHOOL LINEAR ACCELERATOR (LINAC) 

The NPS LINAC is a 100 MeV accelerator technically similar 

to the Stanford University Mark III linear accelerator which 
is described by Chodorow e_t ad. [51] . A detailed description 
of the NPS LINAC is given by Barnett and Cunneen [52] and a 
schematic floor plan is given in Figure III-l. The accelerator 
receives up to 22 megawatts of peak power from each of the 
three klystrons. The machine is pulsed at a rate of 60 pulses/ 
sec, each pulse being about 1 microsecond in duration. 

Electrons are injected at 80 kev and can be accelerated to 
energies between 15 and 100 MeV with an average beam current 
of about 5 microamperes. The beam energy is selected by the 
beam deflection magnets . The magnetic field of these magnets 
is measured by an NMR probe. The beam is focused down to a 
diameter of about 2mm at the target by the quadrupole magnets . 
The spread in the beam energy is controlled by the slit 
settings . 

B. THE SPECTROMETER AND COUNTING SYSTEM 

The- spectrometer is a 120°, 16 inch double-focusing 
magnetic spectrometer described in detail by Oberdier [44] . 

Ten detectors span a portion of the focal plane of the spec- 
trometer, thus giving 10 counting channels for each spectro- 
meter setting. Each detector (a NE102 plastic scintillator) 
is 7/16 inch wide and 1/16 inch thick. The 7/16 inch width 
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allows a momentum spread Ap/p of about 0.3% for each counter. 
Behind these ten detectors are two large detectors, both of 
which span the same portion of the focal plane of the spec- 
trometer as do the 10 smaller scintillators. An electron 
passing through one of the ten scintillators produces a light 
pulse, which is collected by the photomultiplier tube, ampli- 
fied and transmitted to the control room where it is discrimi- 
nated, and is then sent to the scaler if both backing counters 
simultaneously registered outputs. This triple coincidence 
pulse to the scaler must also occur during the time that the 
master gating circuit is pulsing the scaler in order that the 
count will register. The gating circuit prevents background 
counts from registering between beam pulses. A block diagram 
of the counting system is given in Figure III-2. 

The scalers in channels 4 through 7 are capable of re- 
cording periodic pulses up to about 150 M^. while the other 

z 

six channels of the ten channel ladder are connected to scalers 

that can record periodic pulses up to 20 M^. . Only data in 

z 

channels 4 through 7 were analyzed since these could accumulate 
data more rapidly. The count rate in these channels was held 
below 60. counts/second on all runs. Count rate correction 
data were obtained and were used to provide the count rate 
correction factor given in Equation (IV-11). The correction 
factor was linear out to about 100 counts/second, thus main- 
taining the counts below 60 counts/second was more than 
adequate. 



59 




Figure III-2 
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The gas target is covered in detail in Appendix D. Both 
H and d 2 were maintained at liquid nitrogen temperature and at 
about 150 psig in order to have a high density of target 
nuclei . 

C. METHOD OF OBTAINING DATA 

Consider first the data which were obtained during a 
"run" to determine the ratio at a given q^ and scat- 

tering angle. The same data were taken for both deuterium 
and hydrogen and consisted of that which was recorded auto- 
matically and a lesser amount that was recorded manually. 

That which was recorded automatically consisted of any- 
where from 30 to 60 individual units of data, depending on 
the run. Each unit of data was called a "point" and consisted 
of the data accumulated at a given spectrometer setting while 
a given number of incident electrons passed through the gas 
target. In each run, the spectrometer was initially set such 
that the scattered elastic spectrum was well above channel 4. 
The first of the data points was recorded there. Then the 
energy of the spectrometer was lowered between 0.02 to 0.05 
MeV, depending on the- energy range the spectrometer was 
sweeping, and another data point was recorded. This process 
continued until the spectrometer setting was below the peak 
of channel 7. The last data point was then recorded. 

After the selected amount of beam charge had been accumu- 
lated at each data point, the following data were automatically 
fed to a teletype machine that both printed a hard copy and 
punched a paper tape, a) the number of counts recorded in 
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each of the 10 channels, b) the number of counts in coinci- 
dence in the two backing counters, c) the time required to 
obtain the given charge, d) the incident energy of the 
electrons, e) the capacitance of the capacitor used to in- 
tegrate the incident beam, f) the setting of the spectrometer 
in MeV, g) the maximum voltage permitted on the capacitor, 
and h) a number which, when divided by 30 (the full scale 
reading of the integrator) , gave the fractional value of the 
maximum voltage at which the integration of the beam current 
stopped. 

The data which were recorded manually included Secondary 
Emission Monitor (SEM) efficiencies and the gas pressure. 

The SEM efficiencies were taken just prior to the first data 
point, at the peak of each channel and after the last data 
point. The gas pressure was recorded at the start and finish 
of each liquid nitrogen "fill" (as defined in Appendix C) , 
which generally occurred about once at each data point. The 
leakage current in the integration circuits was also recorded 
manually and was used to determine the ground loop errors, but 
was not used to calculate the cross section as it was always 
such a small correction. 

The above data were always taken under conditions that 
were as identical as possible for both the deuteron and 
proton. That is, the experimental geometries, the accele- 
rator beam current, the pressure of the target gases, the 
scattered electron energy, the detection and electronics 
systems and all other possible variables were maintained, as 
nearly as possible, the same. 
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2 -2 

Runs were made at values of q that went from 0.025 fm 
_2 

to 0.40 fm . However, scattering from the target windows or 

_2 

a contaminant in the gas masked the data at 0.025 fm , while 

high background levels and small cross sections made the data 
-2 

taken at 0.40 fm unreliable. Thus the reproducible data 

-2 

reported here covered a region from 0.05 to 0.35 fm 

2 

Within this range of q , forward angle scattering was prefer- 
entially selected to emphasize the electric interaction. 

The most significant change in running conditions between 
the deuteron and proton was the change in beam energy made 
so that the scattered electrons would have the same energy in 
both cases. This insured that the transmission of electrons 
through the spectrometer and the efficiencies of the detectors 

would be the same while taking and H ^ data. This also 
2 

meant that the q for hydrogen was slightly different than 
for deuterium. This difference was accounted for in Equation 

11-42. 

The spectrometer was saturated by increasing the current 
in the coils to a setting corresponding to 100 MeV before data 
were taken in either or in order to insure that the same 
hysteresis loop could be obtained on both spectrums. The 
energy defining slits were generally set to give a full width 
at half maximum to beam energy ratio of about 0.4%. The count 
rate was held below 60 counts per second during all runs. The 
number of counts recorded varied from as much as 170,000 to as 
few as 10,000 in a single channel. The number of counts ob- 
tained depended primarily upon the size of the cross section. 
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An empty target measurement was made at the last spec- 
trometer setting as well as at several others throughout the 
energy range in which data had been taken. Usually those four 

or five readings indicated a "flat" background. Only at the 
2 -2 

lowest q of 0.05 F did any structure show itself in these 

empty target background measurements. Thus with all but the 
2 -2 

q = 0.05 F data runs, the full target background measure- 
ments taken at energies well above the peak of each channel 

were averaged and then projected at those same values through- 

' 2 -2 

out each channel's elastic peak. The q = 0.05 F data were 
handled differently because there was scattering from con- 
taminants in the gases and the windows which complicated the 
background spectrum. In that case empty target background was 
taken at each and every spectrometer setting of both the hydro- 
gen and deuterium. These empty target points were then scaled, 
point by point, to the full target background in order to ob- 
tain the full target background under each elastic peak. 

A significant improvement in procedure over other work 
which had been done at this facility is the procedure that 
steps each detector completely across each elastic peak. 
Formerly, the data from all ten channels were used to construct 
one elastic peak by measuring the dispersion between detectors 
and the efficiency of all others with respect to a chosen re- 
ference detector. This calculation of the efficiencies intro- 
duced errors around 0.5% into the former method which were 
avoided here. 
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IV.. DATA ANALYSIS 



A. GENERAL METHOD OF DATA REDUCTION 

The treatment of the experimental data (see Section III-D) 
is described here. The first step was to change the form of 
the data from punched paper tape to computer card data. This 
was accomplished by a program that required the data in paper 
tape form to be read into an on-line teletype which was con- 
nected to the NPS IBM 360. This program caused the data to 
be punched out on computer cards. Because of a peculiarity 
in the data readout system of the accelerator, decimal points 
were punched as colons. A small program was written which took 
the punched deck mentioned above as input data and replaced all 
colons with decimal points in a second punched deck. This deck 
was then interpreted and checked with the printed teletype out- 
put. Hydrogen and deuterium data were separated and were used 
as input to FORTRAN programs T0P#F0RM where # ran from channel 
4 through 7. The T0P#F0RM program is listed in Appendix D. 
TOP#FORM was written primarily to calculate the cross section 
measured by a particular channel, but it was also used to 
provide count rate corrected chta and to normalize all points 
to a standard incident charge, since the automatic termination 
of the integration at each spectrometer setting did not always 
stop at the same value. T0P#F0RM then prints out these count 
rate corrected, charge normalized counts for each spectrometer 
setting. This part of TOP#FORM was used to plot these so- 
called corrected counts. For that plotting run other input 
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data which are necessary to calculate a cross section, such 
as radiation corrections, pressure, etc. were given dummy 
values. These count rate corrected and normalized counts for 
all four detectors for both the deuterium and the hydrogen 
were then plotted and the position of each peak and the shape 
of the scattered spectrum was determined by visually fitting 
a smooth curve through the points. A typical plot of such 
data is given in Figures IV-la and IV-lb. The difference in 
energy between each peak and the final spectrometer energy at 
which data were taken was then measured because the radiation 
corrections depend upon Ae, and this value was put into the 
FORTRAN program TOPRADCR which is listed in Appendix E. 
Additionally, the peak energy (which is used as the final 
energy of the electron) as well as the pressure of the hydro- 
gen and deuterium are entered as input data to TOPRADCR. 
TOPRADCR has the virial coefficients of hydrogen and deuterium 
and the composition and thicknesses of the gas target windows 
stored as constants. These constants and the run data men- 
tioned above permit TOPRADCR to calculate the density of the 
gas for each channel and the Schwinger, bremsstrahlung and 
Landau radiation corrections. The rationale behind the radi- 
ation correction calculations is given in Section IV-C. The 
average background was calculated for each channel by adding 
all counts recorded prior to the start of the rise in each 
channel and dividing by the number of spectrometer settings 
which recorded these counts. 
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Figure IV-la. Spectrometer Energy (MeV) 
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Figure IV- lb. Spectrometer Energy (MeV) 



With these calculations done, the Schwinger correction, 
the product of the Bethe and the Landau corrections, the gas 
density, the SEM efficiency and the background counts for 
each channel were typed on computer cards and entered with the 
appropriate hydrogen or deuterium deck as the data for TOP#FORM 
which then calculated the elastic cross section and the sta- 
tistical error in that cross section. The method of calcu- 
lating the cross section is given in Section IV-B. The method 
of calculating the statistical error is derived in Section IV-D. 

Four cross section ratios (one for each channel) were 

formed from the cross sections, and then a weighted mean ratio 

was calculated. The electric form factor of the neutron, 

G„ / was then calculated by using this weighted mean ratio, R , 

Jti m 

n 

and the appropriate cross sections and magnetic corrections 
and form factors as shown in Equations IV-14 and IV-15. The 
values of the proton electric form factors and the proton and 
deuteron Mott cross sections and magnetic corrections were 
calculated from computer programs called TOPROTON and TOPDEUTR 
which are listed in Appendix E. Values of D c were from data 
in Appendix C. 



B. THE CALCULATION OF ELASTIC CROSS SECTIONS 
1. The Experimental Cross Section 

The experimental cross sections have been calculated 
from the expression 



a 

Exp 



N K K K 
sc s b I 2^-1 
- - — T cm ster 

N in \ ^ 



(IV-1) 
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where 



N is the number of scattered electrons counted in a 
sc 

particular channel, 

, N. is the number of incident electrons, 

' m 

2 

77 is the number of target atoms per cm , 

A-n. is the so3.id angle defined by the spectrometer 
slits, and 

K ' 3C and K are radiative corrections (see Section 
s b I 



IV-B) , 



is given by 

p t N 






A cos cp ' 



(IV-2 ) 



where 



and 



P 

t 

N 



9 



is the density of the target (gm/cm ) , 
is the thickness of the target (cm) , 
is Avogadro's number, 

is the atomic weight of the target gas. 



is the angle between the normal to a plane solid 
target and incident beam of electrons. In the gas target 
analysis 9 is set equal to zero as the angle cp is not de- 
fined, i.e., rotation of the gas target still presents the 
same number of target atoms per square centimeter. 



N. is determined by 
m 2 

n. = & = -2L 

m e e£ 



(IV-3) 



SEM 



where 



Q = total charge of incident electrons 

-9 

e = electron charge of 1.60210 x 10 coulombs 
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c = capacitance of integrator in farads, 
v = voltage of integrator in volts, 

and 

e = efficiency of SEM. 

Equation IV-1 is appropriately used to calculate 
absolute cross section, but it was used in this experiment 
even though absolute cross sections were not measured. This 
is a correct procedure because, in this experiment, the ratio 
of the cross sections of the deuteron and proton was formed 
and the systematic errors effectively concelled as all data 
were taken under almost the same conditions for both the 
proton and the deuteron. 

N , the total number of scattered electrons, is 
sc 

obtained by integration. Consider Figure (IV-2) , where the 
number of detected electrons in a particular counting channel 
is plotted against the spectrometer energy. 

Let E^ be a spectrometer setting 

and 



the number of counts recorded in this channel 



at that spectrometer setting. 



Then 



N 



sc 



area under scattering peak 

(in Counts - MeV) 

5e 



(IV-4 ) 



where 

6e = the energy width of an individual detector. 

This is because each detector has a finite dimension 
and therefore counts all electrons in the spread Ae. Thus 
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Spectrometer Energy (E) 



Figure IV-2. Spectrum of Counts in a Single Detector 



in obtaining N gc we have: 
E f 



N 



SC 



= 1 



dN 

dE 



dE where 



counts per interval of energy, 
approximated by the sum 



n-l %(N + N ) 

-k-t V' 1 J+I 

N sc = 6 e 



is the number of (IV-5) 
This integral is numerically 



(E i- 



E i+ i) 



(IV-6) 



The energy width 6 e is obtained in terms of E^ by using the 
dispersion relation for the spectrometer optics, 



6e co 

E " Dr 



(IV-7 ) 



o 

where D has been experimentally determined to be 3.92 [53], 
r Q is the mean radius of the spectrometer, 16 inches, and 
1 W is the width of each detector (NE102 plastic scintilla- 
tors) which is 7/32 inches. 
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Therefore 



0.0034877 . 



(IV-8) 



PE _ 7/32" _ 

E ± _ 3.92 (16" ) “ 

Putting these together, we have 

N sc = hi (Ei-E i+1 )/0. 0034877 E.. 

The energy E^ in the denominator is the spectrometer energy 
setting. This setting is very near the center energy of the 
band of electron energies which the spectrometer detects when 
set to the reading E^. Hpwever, the energy of the electrons 
which hit counter 10 will be higher than any other detected 
electrons and the energy of the electrons which hit channel 1 
will be the lowest in energy. In order to take into account 
the fact that each detector is recording electrons with dif- 
ferent energies them the spectrometer setting E^ , several runs 
were made to determine both the absolute calibration of the 
spectrometer and the relative energy spacing of each of the 
four center channels. By comparing the calculated elastic 
peak energy with the spectrometer setting E^, it was deter- 
mined that the peak of channel number 5 was 0.17% higher than 
the energy setting of the spectrometer. Thus channel number 
5 was chosen as the reference energy channel. By successively 
lowering the spectrometer setting from an energy above the 
elastic peak of channel number 4 through channel number 7, the 
dispersion of the magnet was measured. This procedure was 
done at LINAC energies of 17 MeV, 30 MeV, 60 MeV, and 90 MeV. 
The average dispersion of these measurements showed that the 
average energy of the electrons which counter r\ (where 7] is 
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4,5,6, or 7) detected was given by 

E^ - E ± (1.0 + 0.0039 (77-5) ) (IV-9) 

when the spectrometer is set at E^. 

Finally, 

^ = ct 
d sl Exp 

^ _ ^N i +N i+1 J^E i -E i+1 ^y0.0034877E i (1.0+0.0039 (17-5))^^^ 

(^7 ) ( fi - *? a c:i ' )(l.831 x 10 -3 steradians) 

Vee SEM y V A ^ ' 

( iv-io) 

_3 

where the solid angle, A _rv is 1.831 x 10 steradians. 

2 . Corrections to the Observed Counts 

a. Equation IV-9 would correctly measure an absolute 
cross section if: 

i) N., the number of counts recorded on the 
i 

scalars, was the actual number of electrons scattered, i.e., 
if the electronic system did not have a finite resolving 
time. Even though the scalars are capable of counting peri- 
odic pulses up to rates of 150 MHz, they lose counts when 
random (Poisson) pulses are counted at high rates. At 60 
counts per second, which was the highest count rate which was 
permitted during this experiment, about 6 % of the counts were 
lost. 

ii) consisted only of counts due to the 
scattering from the target gas and was not contaminated by 
background counts. 

iii) The efficiency of the detectors were known. 
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b. The count rate correction was determined experi- 
mentally by scattering from an aluminum target at count rates 
up to 500 counts/sec. The number of counts detected was 
plotted as a function of the count rate. Up to about 100 
counts/second the number of counts recorded decreased linearly 
with rate. A linear fit was made to the data in this region, 

^ ^ve N(count rate corrected) = (IV " 11) 

where N/t is the average count rate determined by dividing the 
recorded number of counts by the time of the integration. 

c. The background was determined by counts obtained 
at energies greater than the elastic peak of either H 2 or D 2 / 
whichever was the target at that time, for each of the four 
detectors. The average number of counts for each channel, in 
this region, was then determined and was used as the background. 
During each run an empty target background was also taken to 
insure that there was no rise or fall in this extrapolated 
average background under the elastic peak caused by some source 
other than the or D^. 

d. The efficiency of each detector was not deter- 
mined because the object of this exxueriment was to obtain the 
ratio of the deuteron cross section to the hydrogen cross 
section. Because this ratio was determined for each of the 
four detectors separately, the efficiency of each detector can- 
celled and therefore is not included in the determination of 
the cross section. 

e. An electrometer was used to measure the inte- 
grated incident beam current for each data point. The cir- 
cuitry associated with the electrometer was set to stop the 
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integration and the counting system whenever a predetermined 
voltage across a capacitor charged by the current from a SEM 
was reached. The electrometer did not always stop at exactly 
the same voltage but varied sometimes by as much as a percent 
or two. Therefore both the hydrogen and the deuterium counts 
were normalized to a standard voltage ( v Re f Grence ) f° r eac h 
run. Thus, if v Actua ^ was the actual voltage at the end of 
the integration, then the number of counts recorded was cor- 
rected by 

N (corrected for voltage) = Ih . , . (IV-12) 

c V Actual 

The cross section, after correcting for count rate, 
background and voltage cutoff is 



d£ 
d a 




V 



N. 



ref 



i V 



act 



(l-0.00097Rate i ) 



BKG 



V 



ref 



N 



i+1 



v 



act 



(1-0 . 00097Rate i= _j) 



BKG 



(S-) (- 



SEM 



fi t Na cm 
A 



)(i. 



831 x 10 



-3 



steradians . 



T 

0.0034877 (1.0+0039 (r] 



-5) ) E . )" 



(XV-13) 



C. RADIATION CORRECTION CALCULATIONS 

The proton cross section was derived by Rosenbluth [26] 
with the assumption that a single virtual photon is exchanged 
in the scattering process as shown in Figure IV-3. 

The deuteron cross section, as derived by Jankus [13] in 
the impulse approximation, also assumes that one virtual pho- 
ton is exchanged as shown in Figure IV-4. 
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Feynman Diagram of Rosenbluth Calculation 
Figure IV-3 




Feynman Diagram of Impulse Approximation 
Figure IV-4 

Unfortunately, when one does an electron scattering 
experiment these processes are not the only ones which occur, 
and the experimenter must take into account competing pro- 
cesses in order to be able to compare experiment and theory. 

In fact the experimental scattered energy spectrum which is 
measured contains the effects of three such competing pro- 
cesses which must be accounted for before one can obtain the 
cross sections for processes as shown in Figures IV-3 or IV-4. 
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These competing processes are: 



1) Photon emission and absorption while the electron is 
in the field of the scattering nucleus of interest. The .cor 
rection to the cross section due to this process was first 
calculated by Schwinger [54] and is commonly called the 
Schwinger Correction. 

2) Bremsstrahlung emission in the target before or after 
the nuclear scattering of interest. 

3) Excitation and ionization of material in the target. 

Each of the above processes causes the energy of a 

scattered delectron to be less than it would be if the only 
process which occurred was the one photon exchange. Thus, a 
radiative tail appears on the low energy side of the elastic 
peak, as shown in Figure IV-5. 



E 



f 



d£ 

d-ndE 




Energy of Scattered Electrons 



Elastic Peak with Radiation Tail 



Figure IV-5 
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To obtain a cross section for the nuclear scattering 
event of interest, the experimenter would have to continue 
down the radiative tail to zero energy in order to capture 
those scattered electrons which have some probability of 
losing all their energy in the three processes listed above. 

In practice, the measurement is stopped at Ae below the 
elastic peak (see Figure IV-5) . The value of Ae is generally 
of the order of at least four half-widths (T/2) [55] . Thus 
fewer electrons are detected and the experimental cross 
section is less than it should be. Since each of the three 
processes mentioned above contributes to the radiative tail, 
corrections must be made to account for the undetected elec- 
trons with energy less than - Ae. 

1. Schwinger Correction 

The largest of the corrections is that which accounts 
for photon emission and absorption while the electron is in 
the field of that nucleus which causes the large angle scat- 
tering which is of interest. This correction to the cross 
section was generally about 15% in this experiment. 

The Schwinger correction accounts for many inter- 
actions which may occur in addition to the one virtual photon 
exchange shown in Figures IV-3 and IV-4. Several of these 
possible interactions are shown in Figure IV-6. 

In processes (a) and (b) real photons are emitted 
before and after the nuclear scattering. Processes (c) , (d) , 

and (e) are called vertex renormalization terms and (f) is the 
vacuum polarization term. Improved radiative corrections to 
electron-proton scattering have been calculated by Tsai [56] 
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Feynman Diagrams of Schwinger Correction Interactions 
Figure IV-6 



and by Meister and Yennie [57] and include the effects of 
proton recoil. Numerical evaluation of the Schwinger cor- 
rections in this work used the results given in the review 
article by Maximon [58] * 




where E^ and are the electron energies prior to and after 

scattering, cc is the fine structure constant, m is the elec- 
2 

tron mass, q is the four-momentum transfer, and 



80 



v*> » - r « 

o 

is the Spence function. 

2 . Breinsstrahlung Correction 

When passing through the target material, the elec- 
trons also radiate away energy through radiative collisions 
with atomic electrons and nuclei other than the target nuc- 
leus of interest. The energy loss due to bremsstrahlung has 
been studied by Bethe and Heitler [59] , who calculated the 
cross section $ (E q ) and the average energy loss due to this 
process. By approximating this cross section Heitler [60] and 
Bethe and Ashkin [61] obtain: 




where X Q is the radiation length, which is the distance over 
which the electron has its energy reduced by a factor 1/e. 

In this experiment, the expression given by Butcher and 
Messel [62] has been used to calculate the radiation length 

V 

„ A 

X """* o _ t / o crn 

4 QN q r^z (z+0 . 8) p 0n( 183z 7 ) 

where N q is Avogadro’s number and r^ is the classical elec- 

3 

tron radius, p is the density in gm/cm and A is in grams/mole. 

The target, which consisted of the hydrogen or deu- 
terium gas, was in a pressure vessel which had 0.001 inch steel 
windows for the incident and scattered beam. For the brems- 
strahlung correction, and also the Landau corrections, the 
target was treated as a homogeneous molecular compound of 
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steel and either 11 ^ or which consisted roughly of Fe^ 

/ or Fe.. D~ r , where the “molecular formula" was calcu- 
2 o 1 2b 

lated for each gas pressure. This formula accounts for the 
1 mil windows and the 2 inches of target gas at 150 psig and 
77°K. Following Butcher and Messel the radiation length for 
such a compound target was calculated from 




In the actual FORTRAN programs the X q were multiplied 

2 

by the average density of this compound to give X q in gm/cm . 

Bethe and Heitler [59] also give the probability that 
an electron traversing t radiation lengths will have an energy 
loss less than a times its initial energy (aside from the 
ionization loss and the scattering from the nucleus of 
interest) . That probability is 

a (t/«" 2) 

P (a) - T(l+ t/Qn 2) • 

Hofstadter has integrated this and calculated the cor- 
rection factor to the observed cross section; this brems- 
strahlung correction factor may be written as [6] 

= exp (fift ) log ( 3 e J ' 

o 

where T is the thickness of the target, including the two 

2 

steel 0.001 inch windows as well as the gas, in gm/cm , and 
E is the geometric mean electron energy while going through 
the target. was about 1.02 in this experiment. 
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3 . Ionization Correction 



As the electrons in the beam pass through the target 
they also lose energy by excitation and ionization in the gas 
and in the gas target windows. Landau [63] first calculated 
both the average energy a monoenergetic beam would lose, and 
the spread, or straggling, of the energy distribution of the 
beam after passing through a thin target. From this distri- 
bution it is possible to calculate the fraction of electrons 
which are not detected because data is taken only down to 
E q -Ae. Figure IV-7 shows a theoretical spectrum of monoener- 
getic incident electrons after passing through a layer of 
material. r is the full-width at half -maximum. Qp ro b 
most probable energy loss due to ionization and Ej NC is the 
incident energy, and N(E) is the number of electrons detected 
at energy E. 




Effect of Ionization on Spectrum of 
Monoenergetic Incident Beam 

Figure IV-7 
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Following Landau's [63] method, but using the general 
notation of Breuer [64], one can calculate the fraction, Y (Ae) , 
of the beam one does not detect (because of ionization losses) 
if data taking stops Ae below the peak in Figure IV-7 as 



2 

¥ (Ae) = ^ N (E) dE / ^ 



E . 
me 



N(E)dE . 



o “o 

Landau [63] gave a graph of Y (Ae) but he also derived 
an approximation for ¥(Ae), as have others [65], which has 
been used in this work: 



Y (Ae) = i . 

where the value of to as provided by B6rsch-Supan [66] is 



to = 



X(l- 



X+A 

X+l 



) 



where A = y - 1 = -0.4228..., 

7 is Eulers constant, 

Ae 

X = - 0.05 (is an approximation to the X as defined 

by Landau and Bergstrom [67], 



2 r\ . z . 

§ = 0.154 y V^A 1 T MeV ' and 
T 'i i 

2 

T = target thickness is gm/cm . 

77 ^ is the number of atoms of element (i) per molecule 
of the target which here was considered to be a homogenous 
mixture of the 0.001 steel windows and the gas in the target 
vessel, and z^ and A_^ are the atomic number and atomic weight 
of element (i) . 

The ionization correction is given by 
K. = (1- Y(AE)) _1 
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In this experiment 3C was less than 1.01, i.e., the 
Landau ionization correction was less than one percent. 

D. ERROR ANALYSIS 

1 . Systematic Errors 

Because this experiment measured the ratio of the 
deuteron to the proton cross section, many systematic errors 
in each of these cross sections, such as solid angle, scat- 
tering angle, SEM efficiency, target thickness, and counting 
rate corrections were effectively cancelled when the ratio of 
the two cross sections were formed. This follows because the 
data for both cross sections were taken during one continuous 
LINAC running period with the geometries the same and the beam 
characteristics, excluding the incident energy, as nearly 
identical as possible. 

A systematic error which does not cancel in this 
ratio experiment is that which can arise because of possible 
impurities in the hydrogen and deuterium gases. The supplier 
of the target gases advertized 99.5% pure deuterium and 
99.995% pure hydrogen. However, after closer investigation 
it was discovered that their analytic equipment can only de- 
tect 0.3% or more contamination of D 2 i n H 2 or H 2 D 2* 
the data taken in this experiment used only one bottle of 
deuterium and one of hydrogen, thus all data were taken with 
the same amount of possible contaminant. If there were a 
difference in the contamination of H 2 or this difference 

would have caused a non zero intercept of the plot of G E 

2 ' n 
versus q . For example, a 0.3% contamination of and a 
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0.0% contamination of 11 ^ would have lowered the intercept 
by about 0.0015. 

Whenever empty target background was compared to the 

background above the peak, the deuterium background always 

had the same spectrum distribution in both the empty and 

2 2 

full cases. However, at lower q (particularly q = 0.05) 
the hydrogen spectrum above the hydrogen peak clearly showed 
a contaminant in the gas which, if it were a low atomic 
number contaminant would have been a few percent, but if it 
were a high atomic number contaminant it would have been less 

than a few tenths of a percent. The closeness of the inter- 
2 

cept of G e versus q to zero (we report 0.0000 ± 0.0012 for 
n 

LF model 15) indicates that the H 2 and D2 contaminants were, 
in fact, almost equal. The following figure shows the dif- 
ference between the H 2 and D 2 spectra. The elastic peaks in 
these cases would be at about 1000 counts. 




H 2 "contaminant" peak 
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2 . Random Errors 

The random errors in this ratio experiment, besides 
the statistical (Poisson) counting errors, came from errors 
in determining the density of the gas in the target, vari- 
ations in the SEM efficiency during a run, energy shifts of 
the incident beam, leakage currents in the beam current 
integrators, and very small errors in the radiation cor- 
rections resulting from errors in determining Ae. 

The statistical error for each cross section, due to 
the Poisson distribution of the number of electrons scat- 
tered for a given number of incident electrons, was calcu- 
lated by a subroutine in the same program that calculated 
the cross section and included the statistical errors made 
in determining the background subtraction. To determine the 
fractional error in the cross section caused by these sta- 
tistical fluctuations, the cross section, X, may be written 

as M 

X = Constant • 2 (N v -B) 
k=l K 

where 

N-^ = Total counts (count rate and voltage corrected) 
obtained at the k^ spectrometer setting, 

B = Average background to be subtracted (determined at 
energies above the peak) and assumed to be the same 
for each spectrometer setting, and 
M = The number of settings of the spectrometer. 

M 

Let N = 2 N v be the total number of counts, then: 
k=l K 

X = Constant • (N-M*B) . 
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The variance of the cross section is 



a 2 = (6X) 2 a 2 (|x) 2 : 

x VbN/ N \bB/ ] 



which gives 



= (Constant) 2 (N+M 2 ^) 

X r> 

2 2 
where a B is the variance in the average background and cr^ = n. 

J th 
Let B t = Eb (where = background at 1 setting) and 

J is the number of spectrometer settings above the peak at 

which background was taken. 

Since B 



SB. 

1 



2 1 J 9 l ^ i 

then % = E a B = 4 2 B.= \ B T 
J 1 J J 



B 

J* 



Thus 

and the fractional standard deviation of the cross section is 



a 2 = (Constant) 2 [N+M 2 -■ ] 
X »J 



|n+m 2 b/j |^ 

[N-MB] 2 



which can also be expressed as 
2 



x 

2 



N 



MB(1+ M/J) 



N 



where N q = N-MB and the first term gives the variance due to 
the true counts and the second gives the variance resulting 
from the background. 

The errors in the measured cross sections caused by 
errors in determining the density of the gas in the target 
varied slightly from run to run, but were always less than 
about 0.1%. The pressure dropped while liquid nitrogen was 
being pumped and then rose again and stayed steady until the 
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next nitrogen fill, but the variations were almost identical 
in the hydrogen and the deuterium, and were about 0.3% of the 
measured pressure. In calculations, the average pressure at 
the start of the nitrogen fill or at the completion of the fill 
was the pressure actually used, depending on which was more 
constant. In calculating the gas density, of course, if the 
pressure at the start (completion) of the fill was used for 
hydrogen, then the pressure at the start (completion) of the 
fill was also used for deuterium. 

Generally, six SEM efficiencies were measured during 
a data run on each gas; one was taken at the start of the run, 
one at each of the four peaks, and one at the completion of 
data taking. At each point, four or more measurements were 
made, and the fractional standard deviation of a single mea- 
surement was usually less than 0.2%, giving a fractional 
standard deviation of the mean of the SEM readings of less 
than 0.08% for the six readings. 

Because of problems such as variations in line voltage, 
the beam energy was not constant and thus the scattered energy 
also would change slightly from time to time throughout a run. 
The result of these changes was a shift of counts from one 
channel to another of the counting ladder. This was most 
noticeable when data were being taken along the steep slopes 
of elastic peaks; often a low count in one channel would be 
accompanied by a high count in an adjacent channel. The 
obvious method of reducing these shifts was to close the 
energy defining slits more, but this increased the background 
and thus the ratio of real counts to background became 
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less - giving larger errors on the measured cross sections. 
Nevertheless, the experimenter could not always unambiguously 
estimate how many counts were added or subtracted from a given 
elastic peak due to energy shifts. 

An examination of many runs shows that there was no 
apparent energy shift problem in those runs where the cross 
section was large and the slits could be closed and still 
maintain a good peak to background ratio. In other runs 
where the slits were opened wider, it appeared that a varia- 
tion of as much as a percent could have occurred at several 
points in an elastic peak. Adjacent data points and general 
elastic peak shapes were used to smooth out these large fluc- 
tuations. It is estimated that in the poorer runs these 
energy shifts could have caused an error in the average mea- 
sured cross section ratio of as much as 0.2%. 

There were slight leakage currents in the circuits 
used to measure and integrate the beam current. These cur- 
rents were reduced by improving the insulation of various 
components and by cleaning all sensitive elements. The 
leakage currents on the SEM and the faraday cup, along with 
all the circuits used to accumulate the currents, were mea- 
sured and found to be almost trivial. These measurements 
were repeated during each run and remained very constant. If 
the electron beam was perfectly constant during all parts of 
a given run, the leakage current would introduce no errors. 
Usually the electron current was held constant to within 20%; 
this would cause an error in the ratio of the cross sections 
of about 0.02%. 
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Finally, there is the error in measuring the exact 
position of the peak of the elastic spectrum from the plotted 
data. The amount of this error varies with such things as 
beam stability, energy, and the widths of the slits but 
generally was about 0.01 MeV. This error causes an error 
in the measurement of Ae which in turn causes an error in the 
Schwinger, bremsstrahlung, and ionization corrections. The 
error in the measured cross section due to this Ae error for 
all three radiation corrections was of the order of 0.05%. 

The error in Ae also causes an uncertainty in the final energy 
value at which one terminates the integration of the elastic 
spectrum. This error caused about a 0.03% error in the mea- 
sured cross section. 

Combining all the above errors except the stati stical 
errors and the energy shift errors (since these two vary so 
much from run to run) , it is seen that the total fractional 
error in the measured cross section is about 0.14%. In 
addition to this, the statistical errors in the cross sections 
varied from an average of less than about 0.2% to a maximum of 
about 0.4% for a few runs, while the energy shift errors varied 
from apparently zero through an average of less than 0.1% with 
about 4 runs possibly having errors as large as 0.2%. A sum- 
mary of the errors discussed above is given in Table IV-1. 

The average error in the cross section ratio due to 
all random errors except the statistical error, others’ 

has been calculated from the values of the average cross 
section errors in Table IV-1 to be 0.24%. The total error of 
the weighted mean ratio R was then formed from 
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(iv-1) 






Total 



2 2 

^Stat. 4 S Avg. Others 



TABLE IV-1 

Random Errors in Cross Section Data 
Errors Average Values of Errors 

0.3% 

0 . 1 % 

0.08% 

0 . 10 % 

0.014% 

0.05% 

0.03% 



Statistical 
Gas Density 
SEM Efficiency 
Energy Shifts 
Leakage Currents 

Ae Measurement (Radiation Corrections) 
Ae Measurement (Integration) 



Two significant sources of error were avoided in this 
work that were present in previous work [2]. The use of a gas 
target in place of a polyethylene target and the independent 
determination of the cross section ratios by each channel per- 
mitted these improvements. The gas target scatters electrons 
from the target nuclei all along the beam path, but only those 
nuclei which satisfy the scattering geometry, i.e., those 
which will scatter and recoil through both a collimator which 
is about 6 inches from the gas target and the spectrometer 
slits which are 16 inches from the target, can enter the spec- 
trometer. Thus if the experimenter first scatters from hydro- 
gen, maintains the exact geometry, and then fills the target 
with deuterium and scatters from it, he does not have to 
calculate the target thickness in order to obtain a perfectly 
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accurate ratio of cross sections - if target thickness were 
the only possible error in the cross section determination. 
This is not the case if he used a solid target which is 
thinner than the collimated entrance width to the spectro- 
meter. Stewart [2], who used a thin polyethylene target to 
measure the deuteron and proton cross sections, estimated 
the target thickness error as 0.5%. A second significant 
error which was avoided in this work was the error experi- 
enced in determining the relative efficiencies of each of 
the ten channels. This error was quite informally estab- 
lished to be 0.5% in Ref. 2. In this present work each de- 
tector scanned the complete elastic spectra and the ratio 
of the deuteron cross section to the proton cross section 
was formed separately for the individual detectors. This 
eliminated any detector efficiency errors which would have 
arisen had the method of Ref. 2 been used. 

Having established what errors were made in the 
measured cross sections, one can determine how these errors 



propagated into the value of 



Let the standard deviation 



of a cross section be denoted by and for the proton and 



deuteron. 



Let R = o ‘ d /a p where 



o = experimental cross section of the deuteron and 

oP = experimental cross section of the proton. 

G_, is determined from 
E 

n -V 




1.0 



(IV-14) 
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and G_ /G_ is determined from 
d p 



Cr^ K 

Mott rar 



E 



Mott 



g e 

p 



(IV-15) 



Solving first for the fractional error in R due to 
the errors in both cross sections we have [42] : 



R 

2 



(cr d ) 2 



(aP ) 2 



(IV-16) 



where S„ is the standard deviation of the ratio R. 

R Gth 



The error in the ratio of 



is determined from 



Equation IV-15, thus the standard fractional deviation of the 

o/ 

ratio G /G_ , denoted by S /0 _ /r , , is 

E d E p g e/ g e 

* dp 



j* . h £« 

V% R 



(IV-17) 



From Equation IV-14 we determine the variance in G^ 



to be: 




(IV-18 ) 



and using Equation IV-15, the standard deviation in the 
neutron electric form factor is 

h - '- 2 ' 






1 a Mott K mr 

^ 1 R ^ K P 
Mott m 



w 

D_ 



(IV— 19 ) 



Up to this point in the error analysis only the 



errors in one channel have been discussed, i.e., S_, is the 

R 
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error in R due to errors in the deuteron and proton cross 



sections for any given channel. Letting be the ratio 

measured by channel i and . its standard deviation, then 

Ri 

the weighted mean of all four channels, R , for any given 
run is obtained by calculating 

(IV-20) 




and the standard deviation of this weighted mean ratio is 



S 




(IV— 21) 



R ana are the experimental values of the ratio 

m R^ 

m 

and of the standard deviation of that ratio quoted in 

Table V-I for each run and are the values used to determine 

all derived quantities for each run in this work, e.g., 

2 

G_ , S„ , G, and all others. 

■tj u_, a 



E. REJECTION OF DATA 

Twenty one runs are reported in this work. In each run 
the cross section ratio, R, was measured for each of the four 
channels, giving a total of 84 measurements of R. However, 
nine of these 84 measurements have not been included in the 
data analysis because of clear indicators that each had large 
systematic errors. 
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The data given in Table V-I for q 2 = 0.05 at 0 = 60° 
does not include the data of channel 5 since the " con tami na- 
tion" background subtaction method discussed in Section III-B 
2 

for the q -0.05 data did not work well for that channel. 
Energy shifts caused clear irregularities in four runs. 

For this reason the following were not included in the data 

2 

reported in Table V-I: channel 5 in the second q =0.10 run, 

2 

channel 6 in the first q =0.2 run, channel 7 in the second 
q 2 = 0.2, 0 = 90° run, and channels 4, 5, and 6 in the second 
q 2 = 0.2, 0 = 60° run. 

The inability to establish reliable backgrounds caused 

2 

the rejection of the data of channel 4 in both q =0.35 
runs . 
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V. EXPERIMENTAL RESULTS 



A. NEUTRON CHARGE FORM FACTORS 

The data from 21 runs are reported here. Each run con- 
sisted of four measurements of R, the ratio of the elastic 
electron-deuteron cross section to the elastic electron- 
proton cross section. Of these 84 measurements of R, 9 are 
not included in the reported data for reasons given in 
Section IV-E. 

2 

Table V-I gives the q transferred to the deuteron, the 
scattering angle 9, the incident energies of the electron 
incident on the proton, E. , and deuteron, E. , the experi- 
mental ratio R and the standard deviation of the ratio R as 



determined by the total error (defined in Section IV-D) and 
the statistical (Poisson) error for that run. Also given are 
quantities which were determined by calculations. These 
quantities are explained in Section II. Using the experi- 
mental results and the calculated quantities given in Table 
V-I, the ratio of the deuteron electric form factor to the 
proton electric form factor has been derived and is listed 
in the right hand column. 

The neutron electric form factor G_ was obtained by 

E n 

using Equation 11-43 in the following form 



m 




C 



1 



'E 
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The deuteron charge form factor C E listed in Appendix C and 

2 

the G„ /G and G (q-,) values given in Table V-I were used 
*d E p P • 

to solve for G . The value of C depends upon the deuteron 
n 2 

wave function as well as q . Three sets of the 21 values of 

G_ have been calculated for Lomon-Feshbach models #1, #5, 

E n 

and #15. 

The best estimate of the error in each of these 21 values 

of G e , for each LF model, is the error obtained using s Tota ^ 
n 

from Equation IV-1. The best estimate of the value of 

— ^7 G reported for each of the three models was obtained 
dq 2 E n 

from these sets of G E weighted by s Tota ]_* For comparison, 
d ^ 

7 T G was also calculated from the G„ data weighted by its 

dq 2 E n E n 

statistical error only, and also with equal weighting given 

to each G E to provide an insight into the sensitivity of the 
n 

result to weighting of the data. 

After calculating each G_ and its S_ .. an average G_ 

t iotai E 

2 n n 

and its s To ^- a ^ each q was calculated. A plot of the 

2 

average G E (with its s To -|- a i) versus q is provided for each 
n 

of the three models. For LF model 15 the value of each G^ 

E 

n 

and its S Tota ;L is also given. 

For the 21 rums reported, all values of G £ were calcu- 

n 

lated from the FORTRAN programs listed in Appendix E. In 

addition, in order to check these programs in 15 runs, the 

values of G_ were also calculated on a Wang Model 700 cal- 
E n 

culator using programs not written by the author. These 

Wang, calculations used the same raw data, but used different 

calculations and independent smoothing of the elastic peaks 

to determine G . The values of G determined by each of 
^n L n 
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TABLE V-I (Measured Ratio 
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these methods showed very good agreement, well within the 
statistical errors. This agreement is convincing evidence 
that the computer programs used in the data reduction did not 
contain errors which could bias the conclusions reported in 
this paper. The ratio R reported in Table V-I is the average 
R for those 15 runs. 

2 

The average value of at each q , using Lomon-Feshbach 

n 

Model #1, is given in Table V-II. 

2 

The average value of at each q , using Lomon-Feshbach 

n 

Model #5, is given in Table V-III. 

The value of G calculated from L-F Model 15 and its 
n 

S_ for each of the 21 runs is given in Table V-IV. 

g e 

n 

2 

The average value of G„ at each q , using Lomon-Feshbach 

n 

Model #15 is given in Table V-V. 

B . ROSENBLUTI-I PLOT 

2 

At q w = 0.2 deuteron and hydrogen data were taken at 60, 
75, 90, and 120 degrees in order that a Rosenbluth Plot could 
be constructed. The theoretical interpretation of the 
Rosenbluth Plot is discussed in Section II-G. Table V-VI 
gives the experimental data and the derived quantities used 
as data for the Rosenbluth plot given in Figure V-5. 

C. CHARGE RADIUS OF THE DEUTERON 

The deuteron charge rms radius was calculated by Method 
II given by Schumacher and Bethe [44] . They have provided 
the charge radius of the deuteron from several deuteron wave 
functions. L-F #1 and L-F #5 have a charge radius of 1.971 fm 
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TABLE V-II 



2 

Average G £ at Each q Using Lomon-Feshbach Model #1 
n 



2 

^ - 2 
(fm ) 


C E (LF#1) 


g e 

n 


S„ (Total) 

g e 

n 


0.05 


0.96881 


0.0002 


0.0015 


0.10 


0.93967 


0.0036 


0.0011 


0.20 


0.88667 


0.0028 


0.0011 


0.25 


0.86247 


0.0070 


0.0027 


0.30 


0.83961 


0.0073 


0.0014 


0.35 


0.81796 


0.0095 


0.0019 



A linear fit to the data in Table V-II gives: 



G e = (-0.0003 + 0.0012) + (0.0245 + 0.0059)q 2 , (V-l) 



while a linear fit of G weighted by the statistical 

n 

only gives: n 



G_ = (0.0001 + 0.0010) + (0.0231 + 0.0049)q 2 , (V-2) 



and the 21 points, if equally weighted, give 



G = (-0.0004 + 0.0016) + (0.0263 ± 0.0072)q 2 (V-3) 

E 

n 



A plot of the data in Table V-II is given in Figure V-l. 
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TABLE V-III 



2 

Average G jn at Each q using Lomon-Feshbach Model #5 



q 

(fm“ 2 ) 


C E (LF#5) 


g e 

n 


S„ (Total) 

g e 

n 


0.05 


0.96881 


0.0002 


0.0015 


0.10 


0.93966 


0.0036 


0.0011 


0.20 


0.88665 


0.0028 


0.0011 


0.25 


0.86244 


0.0070 


0.0027 


0.30 


0.83956 


0.0074 


0.0014 


0.35 


0.81790 


0.0096 


0.0019 



A linear fit to the data in Table V-III gives: 

G = (-0.0003 + 0.0012) + (0.0249 + 0.0059)q 2 (V-4) 

E n 

While a linear fit of weighted by the statistical S G 

n E r 

only gives: 

G„ = (0.0001 + 0.0010) + (0.0234 + 0.0049)q 2 (V-5) 

E n 

and the 21 points, if equally weighted, give 

G = (-0.0003 + 0.0016) + (0.0266 + 0.0072)q 2 . (V-6) 

E n 

A plot of the data in Table V-III is given in Figure V-2. 
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TABLE V-IV 



and S (Total) for Each of the 21 Runs 
n & 

n 

Calculated from Lomon-Feshbach Model #15 



q 

(fm -2 ) 


C e (LF#15) 




S_ (Total) 
°E 

n 


0.05 


0.96899 


0.0003 


0.0019 






-0.0005 


0.0023 


0.10 


0.94002 


0.0079 


0.0021 






0.0047 


0.0017 






0.0011 


0.0024 






-0.0012 


0.0021 


0.20 


0.88731 


0.0031 


0.0058 






0.0159 


0.0051 






0.0025 


0.0040 






0.0053 


0.0024 






-0.0002 


0.0024 






0.0023 


0.0025 






-0.0026 


0.0025 






0.0015 


0.0037 


0.25 


0.86324 


0.0061 


0.0031 


0.30 


0.84049 


0.0152 


0.0028 






0.0031 


0.0029 






0.0050 


0.0034 






0.0024 


0.0024 


0.35 


0.81895 


0.0062 


0.0026 






0-.0107 


0.0027 


A plot of 


the data in 


Table V-IV is 


given in Figure 
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TABLE V-V 



Average G^ 

n 

2 


2 

at Each q Using 


Lomon-Feshbach 


Model #15 


Q 

(fm" 2 ) 


C e (LF#15) 


g e 

n 


S„ (Total) 

X 


0.05 


0.96899 


0.0000 


0.0015 


0.10 


0.94002 


0.0034 


0.0010 


0.20 


0.88731 


0.0021 


0.0011 


0.25 


0.86324 


0.0061 


0.0031 


0.30 


0.84049 , 


0.0062 


0.0014 


0.35 


0.81895 


0.0084 


0.0019 



A linear fit to the data in Table V-V gives: 

G = (-0.0000 + 0.0012) + (0.0200 + 0.0058)q 2 (V-7) 

ill 

n 

While a linear fit of weighted by the statistical S„ 

n 

, . n 

only gives: 

G„ = (0.0000 + 0.0010) + (0.0194 + 0.0050)q 2 (V-8) 

E n 

and the 21 points, if equally weighted, give: 

G„ = (-0.0003 + 0.0016) + (0.0231 + 0.0072)q 2 (V-9) 

ill 

n 

A plot of the data in Table V-V is given in Figure V-4. 
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(0.2), Deuteron Form Factor Squared at q = 0.20 fm 
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(0.2) = (0. 7488+0. 0043) + (0. 0013+0. 0013) [1+2 (1+t?) tan^0/2] . (V-ll) 



Qj to 



0.90 



Rosenbluth Plot at q 2 = 0.2 fm 2 
of 

2 

Deuteron Form Factor Squared versus l+2(l+?7)tan (9/2) 



0.85 - 



0.80 * 



0.75 .. 






■^r 



0.70 



0.65 



1+2 (l+T?) tan (0/2) 
Figure V-5 
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while L-F #15 has r 



1.965 fm. The experimental 



~ 2 G E (0) 
dq n 



which were determined using L-F #1, 5, and 15 are plotted 

2 

against the square of the deuteron charge radius, r £ , as 
given by Schumacher and Bethe, in Figure V-6. Schumacher 
notes that such a plot yields an excellent straight line 
for all models with the correct binding energy. The inter- 
section of this line with the value of (0) determined 

dq n 

by the very accurate neutron-electron scattering work yields 
2 

The value of the deuteron charge rms radius was deter- 



E . 



mined from Figure V-6 to be r 



1.9641 + 0.0074 fm. 
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Slope from LF 
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(determined by knowing the wave functions for each model) 



VI. CONCLUSIONS 



From the 75 measurements of the ratio of the elastic 

electron-deuteron cross section to the elastic electron-proton 

2 

cross section at low q the following conclusions are drawn: 

i) the discrepancy in the neutron charge form factor 
measured in several earlier experiments [1,2,3,41 was re- 
solved. The value of ~2 G E (0) was determined to be in 

dq . n 

agreement with the thermal neutron results for all three 
deuteron wave functions (L-F 1,5,15) which were used to analyze 
the data. 

ii) the value of — G (0) ranged from (0.0245 + 

2 dq n 

0.0059) fm when L-F 1 was used to extract G„ from the 

ih 

2 n 

deuteron data, through (0.0249 + 0.0059) fm when L-F 5 was 

2 

used, to (0.0200 + 0.0058) fm when L-F 15 was used. Thus 

L-F 15 produced a value of — G (0) in best agreement with 

*1 n 

the very accurate thermal neutron results of — j G E (0) = 0.0193+ 
2 dc 3 n 

0.0004 fm and was determined to be the model (of the three 

considered) which best describes the deuteron. The method of 
calculating the LF wave functions for different models and 
the more significant differences between the three models con- 
sidered here is discussed in Appendix C. 

iii) the charge radius of the deuteron was measured to 

be r„ = 1.9641 + 0.0074 fm. Prior to this result the most 
E d 

recent measurement [21] was r_ = 1.95 + 0.02 fm. The very 

a 

accurate value of the reported charge radius of the deuteron 
was possible because all data were taken at or below 
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2 -2 

q = 0.35 fm and thus it could all be used to determine the 

2 

intercept and coefficient of the q term in the equation 



G_ (q 2 ) = G (0) - r 2 

E ^ E 6 rms 

n n 

iv) as a check on the results of this experiment a 

2 -2 

Rosenbluth plot of all of the data taken at q = 0.2 fm 

was produced in Figure V-5. The slope of this line is 

(0.0026 + 0.0012). The theoretical slope of this line is 
2 2 

G, ( see Section II-G) . Using the "scaling law" (Equation 

2 M d 2 

11-31) this slope becomes 'J 1 ? ) where G^= (G £ +G E )C E . 

2 P P n 

At q =0.2 the theoretical slope is equal to 0.000822. The 

experimental slope is about 1.5 from this value. 

However, Table V-VI shows, as does Figure V-5, that the 

second 120° point has a significantly different value of G^ 

than do the others. If this one point is dropped, a linear 

fit to the remaining 7 points gives a value of the slope of 

0.0013 + 0.0013 in good agreement with the theoretical slope. 

This shows that the one photon exchange is valid for the 
2 

deuteron at this q . 
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APPENDIX A 



DE VRIES B' FIT 



De Vries et_ al . [10] analysed the proton cross section 

measurements of BumiUer et_ al . [11] which consisted of 58 

data points, the 114 points of Janssens et_ al . [10], and the 

6 points of Chen et_ al , [12] as well as the 71 deuteron cross 

section measurements which de Vries presented. These data 

were used to obtain the best fit to the three pole Clementel- 

Villi [40] model for the isotopic form factors given in 

Equations A-l. This de Vries fit is valid only in the low 
2 

q region. 




15.6 



26.6 





S, 



S 



15.6 



26.6 
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The proton form factors are 



G, 



E. 




P 




while the neutron form factors are 



De Vries started with this theoretical model given by 
Equations A-l and computed cross sections which were then 
compared with the experimental data. The free parameters 

appearing in the model were adjusted so that a minimum in 
2 

X was obtained. 

The mass of the p meson was taken as a free parameter 
since the width of the resonance was so wide (100 MeV) , but 
the seven other free parameters were reduced to six by ap- 
plying the constraint 



Data from the several experiments mentioned above were 
combined in several different combinations. A reduced chi- 
square was calculated for each of these. ‘Combinations b and 
b' were recommended by the authors as the best fits. In this 
experiment the values of the free parameters given by the b' 
fit were used since Buchanan [68] has reanalyzed all known 
proton, neutron, and deuteron data and concluded that 
de Vries b’ fit best describes the proton form factors. The 
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b' combination includes the data of Janssens eit aJL . , de Vries 
et al . , and Chen et; aj. , and gives the following for the 
values of the parameters. 



S„ = 2.628 
E 1 


U) 

!2 

h 

II 


4.193 


S„ =-1.853 
E 2 


cn 

£ 

to 

II 


-3.435 


V = 1.191 
E 1 


I— * 

II 


1.064 


Mp2 = 8.463 







119 



APPENDIX B 



THE VI RIAL COEFFICIENTS OF HYDROGEN AND DEUTERIUM 



In this experiment the target gases were contained in a 
metal gas chamber that was in good thermal contact with a 
liquid nitrogen reservoir which was at atmospheric pressure. 
The chamber was filled from high pressure bottles of H^ or 
D 2 through reducing and regulating valves to a gauge pressure 
of about 150 psig. The density of the gas in the target was 
then computed using the virial coefficients of H 2 and D 2 at 
77.35°K (liquid nitrogen temperature). Since the pressure was 
monitored continuously throughout a data run the virial ex- 
pansion was calculated in terms of pressure coefficients. 

Using the data of Michels and Goudeket [69] and the method 
given in the Handbook of Physics [70] the gas density was de- 
termined to be 



PM 

RT 



1+ B'P + C* 



) gm/ cm^ 



where 

B' = - 0.189409 x 10 -3 in 2 /lb 
C' = 0.1298969 x 10 -7 in 4 /lb 2 . 



3 

At 150 psig hydrogen had a density of 0.0036723 gm/cm 
and deuterium 0.0073391 gm/cm^. 

During a run the pressure was recorded and was then vised 
to calculate (in FORTRAN program TOPRADCR) the gas density 
which in turn was used in Equation IV-2 to obtain the number 
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2 

of target atoms per cm and in the bremsstrahlung and Landau 

radiation correction calculations as well as to partially 

determine the cross section in Equation IV-13. 

The American Institute of Physics Handbook, Third Edition, 

which was published coincidently with the completion of the 

writing of this thesis, has virial coefficients different 

from those given in Reference 70. The data of Reference 70 

was used for all calculations in this work. If the data of 

the AIPH, Third Edition, are used the primary results of this 

work will not be changed. That is, the slope of , the 

n 

radius of the deuteron and the slope of the Rosenbluth plot 

are unchanged. There will, however, be a difference in the 

2 

intercepts of both (0) versus q and of the Rosenbluth 
n 

plot. This change in virial coefficients will cause the 

intercept of G^ to be raised to 0.0023, which might in- 
n 

dicate that the "contaminent" in the discussed in Section 
IV-D is several tenths of a percent. 
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APPENDIX C 



LOMON-FESHBACH DEUTERON ELECTRIC CHARGE STRUCTURE FACTORS 



The deuteron electric charge structure factor, C E , is 
defined as 

C E = $[ u2(r> + v2(r) ] J 0 (f - ) dr 

where u(r) and v(r) are respectively the deuteron S-and D- 

state wave functions normalized so that C^ = 1 when q = 0, 

is the spherical Bessel function of order zero, q is the 

four -momentum transferred and r is the separation between 

the neutron and proton. In order that one might extract the 

neutron form factor from data which give the shape of the 

deuteron one needs to have C„, which means that one must 

E 

have u(r) and v(r). These are determined from the nuclear 
model used. Lomon and Feshbach derived nuclear models from 
nucleon-nucleon data by fitting a boundary condition model 
interaction determined largely by field- theoretic forms [15] . 
Lomon and Feshbach use an energy Independent boundary condition 
for the deuteron wave function. 



F - r fl 
F “ r o br ' 

where r Q is the range of the potential. This condition must 
be different in each angular momentum state [72] . E.L. Lomon 



[14] has provided values of C^ based upon three different 



nuclear models. 



These values of C„ are listed in Table C-I. 

hi 
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The three models which have been considered are called LF#1, 
LF#5 and LF#15 in agreement with Ref. 15. Model L-F #1 has 
a deuteron binding energy of 2.2242 MeV and 4.57 percent 
D-state; model L-F #5 has a binding energy of 2.2245 MeV and 
5.202 percent D-state; model L-F #15 has a binding energy of 
2.2236 MeV and 7.55 percent D-state. 



TABLE C-I 

LOMON-FESHBACH DEUTERON ELECTRIC 
CHARGE STRUCTURE FACTORS 



q 2 (F 2 ) 


c e (lf#i) 


C e (LF#5) 


C E (LF#15) 


0.025 


0.98414 


0.98413 


0.98423 


0.05 


0.96881 


0.96881 


0.96899 


0.10 


0.93967 


0.93966 


0.94002 


0.15 


0.91235 


0.91234 


0.91285 


0.20 


0.88667 


0.88665 


0.88731 


0.25 


0.86247 


0.86244 


0.86324 


0.30 


0.83961 


0.83956 


0.84049 


0.35 


0.81796 


0.81790 


0.81895 


0.40 


0.79741 


0.79735 


0.79851 


0.45 


0.77789 


0.77782 


0.77907 


0.50 


0.75929 


0.75921 


0.76057 


0.55 


0.74156 


0.74147 


0.74291 


0.60 


0.72462 


0.72452 


0.72604 


1.3 


0.54649 


0.54626 


0.54854 


2.0 


0.43328 


0.43293 


0.43555 


3.0 


0.32637 


0.32588 


0.32867 


4.0 


0.25438 


0.25378 


0.25657 
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APPENDIX D 



THE GAS TARGET SYSTEM 

A schematic diagram of the gas target system as it was 
developed for this experiment is shown in Figure D~l. A 
detailed scale drawing of the gas target assembly is given by 
Savage [71] . In this section both the best and the poorest 
techniques which were use'd in this experiment to operate the 
gas target system are discussed and hopefully this will pre- 
vent others from following some unprofitable techniques. 

As the linear accelerator was being turned on for a run 
the gas fill system was connected together between valves 
V7 and V8 and valves VI , V2 and V7 were opened. (Actually 
V7 and V8 were always open during all the runs done with the 
final technique.) The pressure gauge was checked to see that 
it read zero gauge pressure. Then valve VI was closed and 
valves V3, V5 and V8 were opened. The forepump was turned on 
to evacuate the complete gas fill system up to valves V4 and 
V6 which are parts of the gas bottles. This pulled the system 
down to about 50 microns. 

After the accelerator had been tuned for the correct 
energy a large (50 liter) liquid nitrogen dewar was connected 
to feed the liquid nitrogen reservoir. Temperature sensors 
were placed in this reservoir. These sensors were taped to a 
plastic rod and were arranged along this rod as shown in 
Figure D-l. Whenever the nitrogen level fell below the top 
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CD 




Figure D-I 
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GAS TARGET SYSTEM 
SCHEMATIC 




sensor the Cryomiser would open the valve in the compressed 
air line forcing air into the dewar and, since the dewar had 
a pressure seal at the top, this forced liquid nitrogen into 
the gas target reservoir. Through trial and error it was 
found that one minute of pumping was satisfactory; if longer 
"fill" times were set on the Cryomiser it was found that the 
gas target assembly became so cold that it would harden the 
"0" ring seal and a serious vacuum leak would occur. A 
shorter fill time did not adequately fill the reservoir. 

The lower sensor would sound an alarm in the end station if 
the nitrogen level fell below it, and this audio alarm was 
reproduced in the control room. 

In the early runs the gas target was allowed to cool for 
at least 15 minutes before the first gas was put in the target. 
In later runs, the gas was put in immediately after the liquid 
nitrogen was pumping; both methods were satisfactory. 

After a vacuum of at least 50 microns was read on the 
thermocouple and the liquid nitrogen had started to pump, 
valve V2 was closed. A pressure of 50 microns insured that 
any contaminants in the gas fill system would add less than 

5 

one contaminant molecule in 10 target molecules. Both regu- 
lators were turned to the closed position and valves V3 and 
V5 would also be closed. Valve V4 was then opened. The 
pressure gauge marked HP, for high pressure, read between 
2000 to 200 psig depending of course on the amount of gas 
left in the bottle. The LP, low pressure, gauge should then 
read zero. If it did valve V3 was opened. Since the gas 
target had 1 mil stainless steel windows the regulator was 
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moved slowly and steadily in order that these windows would 
not be unnecessarily strained as the pressure increased from 
vacuum to about 150 psig. This pressure was read on both the 
LP gauge and on the much more accurate, larger gauge in 
Figure III-2 simply shown with the word pressure on its face. 
The regulator maintained the system very near the set pressure 
throughout all runs although there usually was some pressure 
drop for a short time while the gas and nitrogen systems 
reached equilibrium. A remote control TV camera v/as panned 
from the gas target to the large pressure gauge after the 
operator had insured that the beam was properly aligned 
through the gas target. The earner a then generally remained 
on the large pressure gauge throughout the run. Since the gas 
pressure would drop whenever a nitrogen fill was progressing, 
the operator recorded the pressure of the gas at the start 
of the fill (i.e., when the audio warning went on in the 
control room) and he would also record the pressure at the 
end of the 1 minute fill (when the audio stopped) . 

After all the hydrogen data were taken valve V4 was 
closed and valves V2 and V5 were opened. Valve V2 was opened 
very slowly in order that the 1 mil stainless window would not 
be unloaded too quickly as it experienced the decrease from 
150 psig. to vacuum. Valve V5 was opened -at this time to 
ensure that no hydrogen leaked through it during the hydrogen 
fill. Again a vacuum of at least 50 microns was pulled on the 
entire system up to valves V4 and V6 (the regulators did not 
provide a vacuum barrier) . After the system was down to at 
least 50 microns, valves V2, V3 and V5 were closed. The 
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regulator on the deuterium was checked to ensure that it 
was closed, thus there was no step-like rise in pressure on 
the gas target windows when valve V6 was opened. Again if 
the deuterium LP gauge read zero, valve V5 was opened and the 
deuterium regulator was turned slowly to reduce the strain on 
the 1 mil steel window. After reading about 150 psig. on 
the large gauge the system had some time to come to equili- 
brium as the operator re-tuned the accelerator. After checking 
the position and appearance of the beam on the target the re- 
mote TV camera was again put on the large pressure gauge and 
the gas pressure was again recorded at the start and finish 
of each nitrogen fill as the deuterium data were taken. 

After taking the deuterium data valve V6 was closed and 
valves V2 and V3 were opened and the system was again pulled 
to at least a 50 micron vacuum. Then the accelerator beam 
was sent through the empty target in order to obtain empty 
target background at several spectrometer settings over the 
range of data taken on that run. 

After the accelerator was secured the forepump and the 
cryomiser were also turned off and valve VI was opened. This 
let the gas system down to air and condensed some gases, e.g. 
water vapor in the gas target. The gas target was not left 
at vacuum because it was believed that the 1 mil windows 
would fail if the target chamber should have had a vacuum 
failure. Very early in these experiments valve V8 was then 
closed in order that excess water would not condense in the 
gas target, but as the liquid nitrogen evaporated during the 
night the 1 mil windows were blown out - presumably due to 
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high vapor pressure within the gas target. The later 
technique was to disassemble the fill line at valve V8. This 
protected the target from a mishap as described above should 
someone have accidently closed V8 or even VI. 

Finally the compressed air line was separated and was 
vented directly into the liquid nitrogen reservoir directly 
above the gas target. This accelerated the rate of evapora- 
tion of the liquid nitrogen and generally would have all the 
nitrogen, and more importantly, all the water which condensed 
during this evaporation removed from the reservoir within 
12 hours. It was found that if the experimenter starts to 
fill the reservoir with liquid nitrogen for another run before 
the condensed water from the last run is completely evaporated, 
then he will most probably cause the indium seal between the 
gas target and the reservoir to fail severely - in fact in 
two experiments the target chamber went to air and only the 
quick closing automatic valve at the collimator saved the 
electron gun. 

It should be mentioned that the gas target with valve V8 
closed will not leak any detectable amount of gas during at 
least a 10 hour period. In the early part of this work valve 
V8 was closed after the initial filling of the target with 
gas, but the pressure at the end of the run generally did 
not reproduce the initial pressure - it was generally higher. 
These pressure variations were of the order of one percent - 
too much of a variation to permit the accuracy of measurement 
required for these experiments. Thus the open target method 



129 



described above was used with continuous monitoring of the 
pressure throughout the run, i.e. by recording the pressure 
at- the start and finish of each nitrogen fill. 
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APPENDIX E 



COMPUTER PROGRAMS USED IN DATA REDUCTION 

The following FORTRAN programs were used to calculate 
quantities explained elsewhere in this work. 

TOPRADCR calculated the Schwinger, bremsstrahlung and 
ionization (Landau) radiation corrections as well as the 
density of the deuterium and hydrogen. 

TOP#FORM calculated the experimental cross section for 
each channel. The sample listed is the program for channel 
number 5 (i.e. 7? =5) . It was also used to obtain the 
smoothed data which was plotted to define the elastic peak. 

TOPROTON calculated the proton's Mott and Rosenbluth 
cross sections of the proton, the magnetic correction term 
and the electric and magnetic form factors of the proton 
based on de Vries' b' fit. 

TOPDEUTR calculated the Mott cross section of the deuteron 
and the magnetic correction term which includes within it the 
relativistic correction to the wave functions. 
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